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ABSTRACT 

Aims. We present abundance analysis based on high resolution spectra of 105 isolated red giant branch (RGB) stars in the Galactic 
Globular Cluster NGC 6121 (M4). Our aim is to study its star population in the context of the multi-population phenomenon recently 
discovered to affect some Globular Clusters. 

Methods. The data have been collected with FLAMES+UVES, the multi-fiber high resolution facility at the ESO/VLT@UT2 tele- 
scope. Analysis was performed under LTE approximation for the following elements: O, Na, Mg, Al, Si, Ca, Ti, Cr, Fe, Ni, Ba, and 
NLTE corrections were applied to those (Na, Mg) strongly affected by departure from LTE. Spectroscopic data were coupled with 
high-precision wide-field UBVI C photometry from WFI@2.2m telescope and infrared 1HK photometry from 2MASS. 
Results. We derived an average [Fe/H] = -1.07 ± 0.01 (internal error), and an a enhancement of [a/Fe] = +0.39 ± 0.05 dex (internal 
error). We confirm the presence of an extended Na-0 anticorrelation, and find two distinct groups of stars with significantly different 
Na and O content. We find no evidence of a Mg-Al anticorrelation. By coupling our results with previous studies on the CN band 
strength, we find that the CN strong stars have higher Na and Al content and are more O depleted than the CN weak ones. The two 
groups of Na-rich, CN-strong and Na-poor, CN-weak stars populate two different regions along the RGB. The Na-rich group defines 
a narrow sequence on the red side of the RGB, while the Na-poor sample populate a bluer, more spread portion of the RGB. In the 
U vs. U - B color magnitude diagram the RGB spread is present from the base of the RGB to the RGB-tip. Apparently, both spec- 
troscopic and photometric results imply the presence of two stellar populations in M4. We briefly discuss the possible origin of these 
populations. 
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1. Introduction 

Observational evidence for variations in the chemical composi- 
tion of light elements in Globular Cluster (GC) stars were known 
since Cohen (1978), who noted a scatter in Na among stars in M3 
and Ml 3. During the last few decades, high resolution spectro- 
scopic studies have definitely confirmed that a GC stellar popu- 
lation is not chemically homogeneous. Even if GC stars are gen- 
erally homogeneous in their Fe-peak element content, they show 
large star-to-star abundance variations in the light elements such 
as C, N, O, Na, Mg, Al, and others (see Gratton et al. 2004 for a 
review). 

During the last twenty years, it has become clear that in red 
giant stars the abundances of some of these elements follow a 
well defined pattern. In particular, there are clear anticorrela- 
tions between the Na and O content, and between Mg and Al. 
Variations in the molecular CH, CN and NH band strengths, 
due to a spread in the abundances of carbon and nitrogen have 
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been observed, as well as anticorrelations between CH and CN 
strengths, and, in some cases, a clear bimodality in the CN con- 
tent. 

Despite the spectroscopic observational evidence collected 
in more than thirty years, the pattern in the light elements is 
not yet well understood. Two scenarios have been proposed to 
explain this observed chemical heterogeneity: the evolutionary 
scenario and the primordial one, both apparently supported by 
observations. 

The observed decline of C content (e.g. in M13, as shown 
by Smith & Briley 2006) and the decreasing ratio 12 C/ 13 C (ob- 
served in M4 and NGC 6528 by Shetrone 2003) along the RGB 
phase, support the evolutionary scenario. According to this the- 
ory, the origin of the observed star-to-star scatter in some el- 
ements is due to the mixing of CNO-cycle-processed material 
transported, in a way not well understood yet, to the stellar sur- 
face. In this way the observed anticorrelations would be present 
in the evolved stages of the life of stars, after the RGB bump. 

At odds with this scenario, in the last few years, spectro- 
scopic studies have revealed light element abundance variations 
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in unevolved main sequence stars and less-evolved RGB stars, 
fainter than the RGB bump. The Na-0 anticorrelation was found 
at the level of the main sequence turn-off (TO) and sub giant 
branch (SGB) in Ml 3 (Cohen & Melendez 2005), NGC 6397 
and NGC 6752 (Carretta et al. 2005; Gratton et al. 2001), 
NGC 6838 (Ramirez & Cohen 2002) and 47 Tuc (Carretta et 
al. 2004). In 47 Tuc, a bimodal distribution in the CN strengths, 
similar to that found among RGB stars (Norris & Freeman 
1979), was found also in the main sequence (MS) by Cannon et 
al. (1998). Moreover, Grundahl et al. (2002) have shown that in 
NGC 6752 the observed scatter in the Stromgren index c\ is due 
to the abundance variations in NH bands in stars both brighter 
and fainter than the RGB bump. This result is consistent with 
a primordial scenario since theory does not predict significant 
mixing below the luminosity of the first dredge-up, observation- 
ally corresponding to the magnitude of the RGB bump. These 
observations suggest that the light element variations should be 
primordial, i.e. they are derived from the chemical composition 
of the primordial site where the GC stars have been generated, or 
alternatively, that a second generation of stars has been formed 
from a medium enriched in some elements (e.g., the self enrich- 
ment model by Ventura et al. 2001). A primordial scenario in 
which such GCs have experienced multiple episodes of star for- 
mation challenges the paradigm that GCs host a single stellar 
population, i.e. that stars of a given cluster are coeval and chem- 
ically homogeneous. 

Very recently, a spectacular, and somehow unexpected con- 
firmation that, at least in some GCs, the origin of the chemical 
anomalies must be primordial came from high precision pho- 
tometry from Hubble Space Telescope observations. The first 
object challenging the paradigm of GCs hosting simple stellar 
populations was to Cen. As shown by Bedin et al. (2004), the 
MS of a> Cen is split into two sequences. But the most exciting 
discovery came from the spectroscopic investigation by Piotto 
et al. (2005), who found that the bluest MS is more metal rich 
than the redder one. The only way to account for the spectro- 
scopic and photometric properties of the two main sequences is 
to assume that the bluest sequence is also strongly He enhanced, 
to an astonishingly high Y=0.38. More recently, Villanova et al. 
(2007) showed that the two main sequences split in at least four 
sub giant branches (SGB) which must be connected in some way 
to the multiplicity of RGBs identified by Lee et al. (1999) and 
Pancino et al. (2000). The metal content of the different SGBs 
measured by Villanova et al. (2007) also implies a large age dif- 
ference among to Cen stellar populations, larger than 1 Gyr. The 
exact age dispersion depends on the assumed abundances (in- 
cluding the He content) of the different stellar populations, and it 
is still controversial (Sollima et al. 2007). Villanova et al. (2007) 
demonstrated that there is also a third MS, running on the red 
side of the two main MSs, and likely connected with the anoma- 
lous RGB-a of Pancino et al. (2000). 

Omega Centauri was well-known since the seventies because 
of its peculiar metallicity distribution. It is the only GC show- 
ing iron-peak element dispersion (Freeman & Rodgers 1975 
and, more recently, Norris et al. 1996, Suntzeff & Kraft 1996). 
In a sense, the findings by Pancino et al. (2000), Bedin et al. 
(2004), Piotto et al. (2005), Villanova et al. (2007) and Sollima 
et al. (2007) could simply be considered additional evidence that 
a) Cen is so peculiar that it might not be a GC. Perhaps, as sug- 
gested by many authors (Freeman 1993, Hughes & Wallerstein 
2000), it might simply be the nucleus of a much larger system, 
likely disrupted by the tidal field of our Galaxy. 
In this sense the most recent discovery by Piotto et al. (2007) that 
the MS of NGC 2808 is split into three, distinct sequences came 



as a sort of surprise, shaking at its foundation our understanding 
of GC stellar populations. NGC 2808 has been always consid- 
ered a GC, with many peculiar properties regarding its metal 
content and its color-magnitude diagram, but a genuine, massive 
GC. Still, it hosts multiple stellar populations. Moreover, also in 
this case, in view of the negligible dispersion in iron-peak ele- 
ments in NGC 2808, the only way so far available to reproduce 
the three MSs is to assume that there are three populations, char- 
acterized by three different helium contents, up to an (again) as- 
tonishingly high Y=0.40. Interestingly enough, D'Antona et al. 
(2005) already made the hypothesis of three groups of stars, with 
three different helium abundances in order to explain the pecu- 
liar, multi-modal Horizontal Branch (HB) of NGC 2808 (Sosin 
et al. 1997, Bedin et al. 2000). Piotto et al. (2007) simply found 
them in the form of a MS split. Piotto et al. (2007) also noticed 
that the different stellar populations in NGC 2808 are consis- 
tent with the spectroscopic observations by Carretta et al. (2006), 
who identified three groups of stars with different Oxygen abun- 
dances. The fraction of stars in the three abundance groups is in 
rough agreement with the fraction of stars in the three MSs. 

NGC 2808 and to Cen are, at the moment, the most extreme 
examples of a rather complex observational scenario. In fact, ev- 
idence of multiple populations has been found in other GCs, like 
NGC 1851 (Milone et al. 2008), NGC 6388 (Siegel et al. 2007, 
Piotto et al. 2008), and M54 (Piotto et al., in preparation) in the 
form of a split in the SGB. In NGC 1851 the presence of a group 
of RGB stars with enhanced Sr and Ba and strong CN bands, 
among the majority of CN-normal RGB stars (Yong & Grundahl 
2008), and the presence of a bimodal HB, agrees with the hy- 
pothesis of two stellar generations inferred by the observed SGB 
split. 

The extremely peculiar HB of NGC 6388 (Rich et al. 1997), 
with the presence of extremely hot HB stars (Busso et al. 2007) 
was already interpreted in terms of multiple, helium enhanced, 
population by Sweigart & Catelan (1998), and in the more de- 
tailed analysis by Caloi & D'Antona (2007). 

As for M54, it has been recognized as the nucleus of the 
Sagittarius dwarf galaxy (Da Costa & Armandroff 1995; Bassino 
& Muzzio 1995; Layden & Sarajedini 2000), and it could simply 
represent what to Cen was long time ago. 

It is worth noting that the clusters showing multiplicities in 
their CMDs are among the most massive GCs of our Galaxy 
(M > 10 6 M o ) and all of them have peculiar HBs, as well as 
peculiar abundances, including Na-O anticorrelations. However, 
we must also note that much less massive GCs, with no evidence 
of multiple populations identified so far, show very large star-to- 
star abundance variations. It is noteworthy to recognize here that 
the Na-0 anticorrelation has been found in about 20 GCs (see 
Carretta et al. 2006 for the most updated list). 

In summary, it is clear that GCs are not as simple systems 
as thought in the past. Up to now, we lack a complete explana- 
tion for the mechanisms necessary to understand the observa- 
tional scenario. A systematic study of the chemical abundances 
of many stars in GCs is needed in order to better understand the 
star formation history of these objects. 

In this work we present a study on the chemical abundances 
of the GC NGC 6121 (M4) from high resolution spectra of its 
RGB stars. 

As far as we know, M4 shows no evidence of multiple stellar 
populations in its CMD, and its mass (log ^ = 4.8, Mandushev 
et al. 1991) is much smaller than the mass of the clusters with 
the photometric peculiarities discussed above. 

Chemical abundances from high resolution spectra of M4 
RGB stars have been already measured by different groups 
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of investigators: Gratton, Quarta & Ortolani 1986 (hereafter 
GQO86), Brown & Wallerstein 1992 (BW92), Drake et al. 1992 
(D92), Ivans et al. 1999 (199), and Smith et al. (2005). These 
authors have found a range of [Fe/H] between -1.3 dex (BW92) 
and - 1 .05 dex (D92). For further details, see 199 who have sum- 
marized the chemical abundances found in previous studies. A 
study by Norris 1981 (hereafter N81) of 45 RGB stars showed a 
CN bimodality in M4, e. g., stars of very similar magnitudes and 
colors have a bimodal distribution of CN band strengths. By an- 
alyzing 4 RGB stars, two selected from the CN-weak group and 
two from the CN-strong one, D92 (and then Drake et al. 1994) 
found differences in the Na and Al content. 199 found that O is 
anticorrelated with N, whereas Na and Al abundances are larger 
in CN-strong stars. Looking at the evolutionary states of these 
stars, both 199 and Smith & Briley 2005 (SB05) did not find any 
strong correlation between CN band strength and the position on 
the CMD. More recently, Smith et al. (2005), studying the fluo- 
rine abundance in seven RGB stars in M4, found a large variation 
in 19 F, which is anticorrelated with the Na and Al abundances. 
In these previous studies, both the evolutionary scenario and a 
primordial one have been taken into account in order to explain 
the light element variations and the CN bimodality in M4. 

In this work, we analyze high resolution spectra in order 
to study chemical abundances for a large sample of M4 RGB 
stars and compare our results with those of previous studies. In 
Section 2 we provide an overview of the observations and target 
sample. The membership criterion used to separate the proba- 
ble cluster stars is described in Section 3, and the procedure to 
derive the chemical abundances is in Section 4. We present our 
results in Section 5, and discuss them in Section 6 and Section 
7. Section 8 summarizes the results of this work. 



m 




B-I 



Fig. 1. Distribution of the UVES target stars on the B vs. (B - I) 
CMD corrected for differential reddening. 



2. Observations and Data Reduction 

Our dataset consists of UVES spectra collected in July- 
September 2006, within a project devoted to the detection of 
spectroscopic binaries in the GC M4, and to the measurement 
of the geometric distance of this cluster (Programs 072.D-0742 
and 077.D-0182). Data are based on single 1200-1800 sec expo- 
sures obtained with FLAMES/VLT @ UT2 (Pasquini et al. 2002) 
under photometric conditions and a typical seeing of 0.8-1.2 arc- 
sec. The 8 fibers feeding the UVES spectrograph were centered 
on 115 isolated stars (no neighbours within a radius of 1.2 arc- 
sec brighter than V + 2.5, where V is the magnitude of the target 
star) from ~ 1 mag below the HB to the tip of the RGB of M4, 
in the magnitude range 10.5 < V < 14.0. 

The UVES spectrograph was used in the RED 580 setting. 
The spectra have a spectral coverage of ~2000 A with the cen- 
tral wavelength at 5800 A. The typical signal to noise ratio is 
S/N ~ 100- 120. 

Data were reduced using UVES pipelines (Ballester et al. 
2000), including bias subtraction, flat-field correction, wave- 
length calibration, sky subtraction and spectral rectification. 
Stars were carefully selected from high quality photometric 
UBVIc observations by Momany et al. (2003) obtained with the 
Wide Field Imager (WFI) camera at the ESO/MPI 2.2m tele- 
scope (total field of view 34 x 33 arcmin 2 ) coupled with infrared 
JHK 2MASS photometry (Skrutskie et al. 2006). Given the spa- 
tially variable interstellar reddening across the cluster (Cudworth 
& Rees 1990; Drake et al. 1994; Kemp et al. 1993; Lyons et al. 
1995; and 199) we have corrected our CMDs for this effect (as 
done in Sarajedini et al. 2007). Fig. Q] shows the position of the 
target stars on the corrected B vs (B-I) CMD. 

3. Radial velocities and membership 

In the present work, radial velocities were used as the 
membership criterion since the cluster stars all have similar 
motion with respect to the observer. The radial velocities 
(v r ) of the stars were measured using the IRAF FXCOR 
task, which cross-correlates the object spectrum with a 
template. As a template, we used a synthetic spectrum ob- 
tained through the spectral synthesis code SPECTRUM (see 
http://www.phys.appstate.edu/spectrum/spectrum.html for 
more details), using a Kurucz model atmosphere with roughly 
the mean atmospheric parameters of our stars T e ff = 4500 K, 
log(g) = 2.0, v t = 1.4 km/s, [Fe/H] = -1.10. At the end, each 
radial velocity was corrected to the heliocentric system. We 
calculated a first approximation mean velocity and the r.m.s (cr) 
of the velocity distribution. Stars showing v r out of more than 
3cr from the mean value were considered probable field objects 
and rejected, leaving us with 105 UVES spectra of probable 
members. After this procedure, we obtained a new mean radial 
velocity of 70.6 + 1 km/s from all the selected spectra, which 
agrees well with the values in the literature (Harris 1996, 
Peterson et al. 1995). 

4. Abundance analysis 

The chemical abundances for all elements, with the exception of 
Oxygen, were obtained from the equivalent widths (EWs) of the 
spectral lines. An accurate measurement of EWs first requires a 
good determination of the continuum level. Our relatively metal- 
poor stars, combined with our high S/N spectra, allowed us to 
proceed in the following way. First, for each line, we selected a 
region of 20 A centered on the line itself (this value is a good 
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compromise between having enough points, i. e. a good statistic, 
and avoiding a too large region where the spectrum can be not 
flat). Then we built the histogram of the distribution of the flux 
where the peak is a rough estimation of the continuum. We re- 
fined this determination by fitting a parabolic curve to the peak 
and using the vertex as our continuum estimation. Finally, the 
continuum determination was revised by eye and corrected by 
hand if a clear discrepancy with the spectrum was found. Then, 
using the continuum value previously obtained, we fit a gaussian 
curve to each spectral line and obtained the EW from integration. 
We rejected lines if affected by bad continuum determination, 
by non-gaussian shape, if their central wavelength did not agree 
with that expected from our linelist, or if the lines were too broad 
or too narrow with respect to the mean FWHM. We verified that 
the gaussian shape was a good approximation for our spectral 
lines, so no lorenzian correction was applied. The typical error 
for these measurements is 3.6 mA, as obtained from the compar- 
ison of the EWs measured on stars having similar atmospheric 
parameters (see Sect. 5.1 for further details). 

4.1. Initial atmospheric parameters 

Initial estimates of the atmospheric parameters were derived 
from WFI BVI photometry. First-guess effective temperatures 
(T e ff) for each star were derived from the T e ff-color relations 
(Alonso et al. 1999). (B - V) and (V - 1) colors were de- 
reddened using a reddening of E(B - V) = 0.36 (Harris 1996). 
Surface gravities (log(g)) were obtained from the canonical 
equation: 

log(g/g ) = log(M/M ) + 4 • log(T eff /T ) - log(L/L ) 

where the mass M/M was derived from the spectral type (de- 
rived from T e ff) and the luminosity class of stars (in this case 
we have a III luminosity class) through the grid of Straizys & 
Kuriliene (1981). The luminosity L/L was obtained from the 
absolute magnitude My, assuming an apparent distance modu- 
lus of (m - M)y = 12.83 (Harris 1996). The bolometric correc- 
tion (BC) was derived by adopting the relation BC - T e ff from 
Alonso etal. (1999). 

Finally, microturbolence velocity (v t ) was obtained from the re- 
lation (Gratton et al. 1996): 

v t (km/s) = 2.22 - 0.322 • log(g) 

These atmospheric parameters are only initial guesses and are 
adjusted as explained in the following Section. 

4.2. Chemical abundances 

The Local Thermodynamic Equilibrium (LTE) program MOOG 
(freely distributed by C. Sneden, University of Texas at Austin) 
was used to determine the metal abundances. 

The linelists for the chemical analysis were obtained from 
the VALD database (Kupka et al. 1999) and calibrated us- 
ing the Solar-inverse technique. For this purpose we used the 
high resolution, high S/N Solar spectrum obtained at NOAO 
(National Optical Astronomy Observatory, Kurucz et al. 1984). 
We used the model atmosphere interpolated from the Kurucz 
(1992) grid using the canonical atmospheric parameters for 
the Sun: T eff = 5777 K, log(g) = 4.44, v t = 0.80 km/s and 
[Fe/H] = 0.00. 

The EWs for the reference Solar spectrum were obtained in 
the same way as the observed spectra, with the exception of the 



strongest lines, where a Voigt profile integration was used. In the 
calibration procedure, we adjusted the value of the line strength 
log(gf) of each spectral line in order to report the abundances ob- 
tained from all the lines of the same element to the mean value. 
The chemical abundances obtained for the Sun are reported in 
Tab. [T] The derived Na and Mg abundances were corrected for 
the effects of departures from the LTE assumption, using the pre- 
scriptions by Gratton et al. (1999). 

Table 1. Measured Solar abundances 

(loge(X) = log(N x /N H ) + 12). 



Element 


UVES 


lines 


OI 


8.83 


1 


NaI NLTE 


6.32 


4 


MgI NLTE 


7.55 


3 


All 


6.43 


2 


Sil 


7.61 


12 


Cal 


6.39 


16 


Til 


4.94 


33 


Till 


4.96 


12 


CrI 


5.67 


32 


Fel 


7.48 


145 


Fell 


7.51 


14 


Nil 


6.26 


47 


Ball 


2.45 


2 



With the calibrated linelist, we can obtain refined atmo- 
spheric parameters and abundances for our targets. Firstly model 
atmospheres were interpolated from the grid of Kurucz mod- 
els by using the values of T e ft, log(g), and v t determined as ex- 
plained in the previous section. Then, during the abundance anal- 
ysis, T e ff, v t and log(g) were adjusted in order to remove trends in 
Excitation Potential (E.P.) and equivalent widths vs. abundance 
respectively, and to satisfy the ionization equilibrium. T e ff was 
optimized by removing any trend in the relation between abun- 
dances obtained from the Fel lines with respect to the E.P. The 
optimization of log(g) was done in order to satisfy the ionization 
equilibrium of species ionized differently: we have used Fel and 
Fell lines for this purpose. We changed the value of log(g) until 
the following relation was satisfied: 

loge(FeII) - loge(FeI) = loge(Fell)* - loge(Fel), 

We optimized v t by removing any trend in the relation abundance 
vs. EWs of the spectral lines. This iterative procedure allowed 
us to optimize the values of atmospheric parameters on the ba- 
sis of spectral data, independent of colors. This is an important 
advantage for those clusters, as M4, that are projected in regions 
characterized by a differential reddening. 

The adopted atmospheric parameters, together with the coordi- 
nates, the U,B,V,I C , and the 2MASS J,H,K magnitudes, are listed 
in Tab. [8] All the reported magnitudes are not corrected for dif- 
ferential reddening. 

Having T e ff determinations independent of colors, we can verify 
whether the T e ff-scale is affected by systematic errors. For this 
purpose we used the T e ff and [Fe/H] of our stars (see Sec. 6.1) 
to obtain intrinsic B - V colors from Alonso's relations. These 
colors were compared with WFI photometry corrected for dif- 
ferential reddening. We obtained a reddening of: 

E(B - V) = +0.34 ± 0.02 

This value agrees well with the E(B - V) = 0.36 of Harris 
(1996) and with 199 who find E(B - V) = 0.33 ± 0.01. 
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Therefore, we conclude that our T e ff-scale is not affected by 
strong systematic errors. 

Finally we present a v t - log(g) relation obtained from our data: 

v t = -(0.254 + 0.016) • log(g) + (1.930 + 0.035) 

It gives a mean microturbolence velocity 0.15 km/s lower than 
that given by Gratton et al. (1996), but this is not a surprise be- 
cause also in other papers (Preston et al. 2006) Gratton et al. 
(1996) was found to overestimate the microturbolence, especial 
in our T e jf regime. We underline that our formula is valid for ob- 
jects in the same T e jf-log(g) regime as out targets, i.e. cold giant 
stars. 



i 1 1 I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 




t, O Eh 



L i i i I I i i i i I i i i i I i i i i I i i i i -I 

6297 6298 6299 6300 6301 6302 

Wavelength [A] 



4.3. The O synthesis 

Instead of using the EW method, we measured the O content 
of our stars by spectral synthesis. This method is necessary be- 
cause of the blending of the target O line at 6300 A with other 
spectral lines (mainly the Ni transition at 6300.34 A). For this 
purpose, we used a linelist from the VALD archive calibrated on 
the NOAO Solar spectrum as done before for the other elements, 
but in this case we changed the log(gf) parameter of our spec- 
tral lines in order to obtain a good match between our synthetic 
spectrum and the observed one. 

With the calibrated linelist, it is possible to establish the O 
content of our targets. To this aim, we used the standard MOOG 
running option synth that computes a set of trial synthetic spec- 
tra and matches these to the observed spectrum. The synthetic 
spectra were obtained by using the atmospheric parameters de- 
rived in the previous section. In this way, the Oxygen abundances 
were deduced by minimization of the observed-computed spec- 
trum difference. An example of a spectral synthesis is plotted 
in Fig.|2]where the spectrum of the observed star #34240 (thick 
line) is compared with synthetic spectra (thin lines) computed 
for O abundances [0/Fe]= +0.30, +0.47, and +0.60 dex. In the 
plotted spectral range, the observed spectrum is contaminated by 
two telluric features (not present in the models), but none affect- 
ing the O line. 

The O line is faint, and we could measure the O content only 
for a sub-sample of 93 stars. In the remaining 12 spectra, the bad 
quality of the Oxygen line prevented us from measuring accurate 
O abundance. 



5. Results 

The wide spectral range of the UVES data allowed us to de- 
rive the chemical abundances of several elements. The mean 
abundances for each element are listed in Tab. [2] together with 
the error of the mean, the rms scatter and the number of mea- 
sured stars (N stars ). The rms scatter (hereafter cr t, s ) is assumed 
to be the 68.27th percentile of the distribution of the measures 
of the single stars and the error of the mean is the rms divided 
by VN stars - 1- In the last column the abundances derived by 199 
were listed as a comparison. 

The derived Na and Mg abundances were corrected for the 
effects of departures from the LTE assumption (NLTE correc- 
tion) using the prescriptions by Gratton et al. (1999), as done for 
the Sun. In the paper all the used Na and Mg abundances are 
NLTE corrected, also were not explicitely indicated. The mean 
NLTE corrections obtained for [Na/Fe] and [Mg/Fe] are -0.02 
dex and +0.07 dex, respectively. For the Fe abundance, we do 



Fig. 2. Spectrum of the star #34240, compared with synthetic 
spectra in the region 6297-6302 A, which includes the Oxygen 
line at 6300.31 A. Synthetic spectra were computed for O abun- 
dances [O/Fe] = +0.30, +0.47, and +0.60 dex. Thick line is the 
observed spectrum, thin lines are the synthetic ones. 
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Fig. 3. Box plot of the M4 giant star element abundances from 
UVES data. For each box the central horizontal line is the me- 
dian of the data, while the upper and lower lines indicate the 
higher and lower lcr value respectively. The points represent the 
individual measurements. 



not distinguish between the results obtained from I and II ion- 
ization stages, because their values are necessarily the same in 
the chosen method. 

Chemical abundances for the single stars are listed in Tab. [9] 
A plot of our measured abundances is shown in Fig. [3] where, 
for each box, the central horizontal line is the median values 
for the element, and the upper and lower lines indicate the lcr 
values higher and lower than the median values, respectively. 
The points represent individual measurements. 



5.1. Internal errors associated to the chemical abundances 

The measured abundances of of every element vary from star to 
star as a consequence of both measurement errors and intrinsic 
star to star abundance variations. In this section, our final goal is 
to search for evidence of intrinsic abundance dispersion in each 
element by comparing the observed dispersion cr obs and that pro- 
duced by internal errors (o" tot ). Clearly, this requires an accurate 
analysis of all the internal sources of measurement errors. We 
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Table 2. The average abundances of M4 stars The results by 199 
are shown for comparison in Column 5. 





This work 






199 










O" obs 


N stars 






[OI/Fe] 


+0.39 


±0.01 


0.09 


93 


+0.25 ± 


0.03 


[NaI/Fe] NLTE 


+0.27 


±0.02 


0.17 


105 


+0.22 ± 


0.05 


[MgI/Fe] NLT E 


+0.50 


±0.01 


0.06 


105 


+0.44 ± 


0.02 


[All/Fe] 


+0.54 


±0.01 


0.11 


87 


+0.64 ± 


0.03 


[Sil/Fe] 


+0.48 


±0.01 


0.05 


105 


+0.55 ± 


0.02 


[Cal/Fe] 


+0.28 


±0.01 


0.04 


105 


+0.26 ± 


0.02 


[Til/Fe] 


+0.29 


±0.01 


0.05 


105 


+0.30 ± 


0.01 


[Till/Fe] 


+0.35 


±0.01 


0.06 


105 


+0.30 ± 


0.01 


[CrI/Fe] 


-0.04 


±0.01 


0.05 


105 






[Fel/H] 


-1.07 


±0.01 


0.05 


105 


—1.18 ± 


0.00 


[Nil/Fe] 


+0.02 


±0.01 


0.03 


105 


+0.05 ± 


0.01 


[Ball/Fe] 


+0.41 


±0.01 


0.09 


103 


+0.60 ± 


0.02 



remark here that we are interested in star-to-star intrinsic abun- 
dance variation, i.e. we want to measure the internal intrinsic 
abundance spread of our sample of stars. For this reason, we are 
not interested in external sources of error which are systematic 
and do not affect relative abundances. 

It must be noted that two sources of errors mainly contribute to 
cr tot . They are: 

- the errors cr EW due to the uncertainties in the EWs measures; 

- the uncertainty <x atm introduced by errors in the atmospheric 
parameters adopted to compute the chemical abundances. 

In order to derive an estimate of cr EW , we consider that EWs of 
spectral lines in stars with the same atmospherical parameters 
and abundances are expected to be equal; any difference in the 
observed EWs can be attributed to measurements errors. Hence, 
in order to estimate cr EW , we applied the following procedure: 

- Select two stars (#33683 and #33946) from our sample char- 
acterized by exactly the same T e ff, and log(g), v t , [Fe/H] 
within a range of 0.05; 

- Compare the EWs of the iron spectral lines in these two stars 
and calculate the standard deviation of the distribution of the 
differences between the EWs divided by v2. We found a 
value of 3.6 mA that we assume to represent our estimate for 
the errors in the EWs. Iron lines were selected because Fe 
abundance has always been found to be the same for stars 
in GCs, with the exception of a> Cen. This means that any 
differences in the EW for the same couple of lines is due 
only to measurement errors; 

- Calculate the corresponding error in abundance measure- 
ments (<r EW ) by using a star (#21728) at intermediate tem- 
perature, assumed to be representative of the entire sample. 
To this aim, we selected two spectral lines for each element 
in order to cover the whole E.P. range, changed the EWs 
by 3.6 mA and calculated the corresponding mean chemical 
abundance variation. This number divided by VNii nes - 1 is 
our o"bw- Results are listed in Tab. [3] 

A more elaborate analysis is required to determine o~ a t m . We fol- 
lowed two different approaches. 

In order to better understand the method, it is important to 
summarize the role of the atmospheric parameters in the deter- 
mination of chemical abundances. As fully described in sections 
4.1, the best estimate of chemical abundances, T e ff, v t and log(g) 
obtained from the spectrum of a single star are those that satisfy 
at the same time three conditions: 



- removing any trend from the straight line that best fits the 
abundances vs. E.P.; 

- removing any trend from the straight line that best fits the 
abundances as a function of equivalent widths; 

- satisfy the ionization equilibrium through the condition 
[Fel/H] = [FeII/H]. 

Therefore, any error in the slope of the best fitting lines and in 
Fel/Fell determinations produces an error on the measured abun- 
dance. 

In order to derive the error in temperature we applied the 
following procedure. First we calculated, for each star, the er- 
rors associated with the slopes of the best least squares fit in 
the relations between abundance vs. E.P. The average of the er- 
rors corresponds to the typical error on the slope. Then, we se- 
lected three stars representative of the entire sample (#29545, 
#21728, and #34006) with high, intermediate, and low T e ff, re- 
spectively. For each of them, we fixed the other parameters and 
varied the temperature until the slope of the line that best fits 
the relation between abundances and E.P. became equal to the 
respective mean error. This difference in temperature represents 
an estimate of the error in temperature itself. The value we found 
is AT eff = 40 K. 

The same procedure was applied for v t , but using the relation 
between abundance and EWs. We obtained a mean error of Av t = 
0.06 km/s. 

Since log(g) has been obtained by imposing the condition 
[Fel/H] = [FeII/H], and the me asures of [Fel /H] and [FeII/H] 
have averaged uncertainties of cr star [FeII/H] and cr star [Fel/H] 
(where <x star [Fe/H] is the dispersion of the iron abundances de- 
rived by the various spectral lines in each spectrum and given by 
MOOG, divided by VNiines - 1), in order to associate an error to 
the measures of gravity we have varied the gravity of the three 
representative stars such that the relation: 



[Fel/H] - cr star [FeI/H] = [FeII/H] + cr star [FeII/H] 

was satisfied. The obtained mean error is A[log(g)]=0.12. 

Once the internal errors associated to the atmospheric pa- 
rameters were calculated, we re-derived the abundances of the 
three reference stars by using the following combination of at- 
mospheric parameters: 

- (T eff + A(T eff ), log(g), V t ) 

- (T e ff, log(g) + A(log(g)), V t ) 

- (Teff, log(g), V t + A(Vt)) 

where T e ff , log(g), v t are the measures determined in Section 4.2. 

The resulting errors in the chemical abundances due to un- 
certainties in each atmospheric parameter are listed in Tab. [3] 
(columns 2, 3 and 4). The values of cr atm (given by the squared 
sum of the uncertainties introduced by each single parameter) 
listed in column 5 are our final estimates of the error introduced 
by the uncertainties of all the atmospheric parameters on the 
chemical abundance measurements. 

We also used a second approach to estimate cr atm as confir- 
mation of the first method. In order to verify our derived un- 
certainties related to the temperature, gravity and microturbo- 
lence, we considered that RGB stars with the same I magnitude 
must have the same atmospheric parameters. Hence, we started 
by plotting the magnitude I as a function of T e ff, log(g) and v t 
as shown in Fig. [4] We see in this figure that the typical inter- 
nal error in our I photometry (~ 0.01 mag) translates into an 
error in temperature of ~ 3 K, in log(g) of ~ 0.005, and in v t 
of ~ 0.002 km/s, each of which is absolutely negligible. The 
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data were fitted by a parabolic curve in the case of T e jf, and by 
straight lines in the cases of log(g) and v t . We determined the 
differences between the T e ff, log(g) and v t of each star and the 
corresponding value on the fitting curve. Assuming the 68.27th 
percentile of the absolute values of these differences as an es- 
timate of the dispersion of the points around the fitting curve, 
all the stars with a distance from the curve larger than 3<x were 
rejected (the open triangles in Fig. The four probable AGB 
stars (filled triangles) that are evident on the B - (B - I) CMD of 
Fig.Q]were also rejected. With the remaining stars, the curve was 
refitted and the differences between the atmospheric parameters 
of the stars and the fit were redetermined. The 68.27th percentile 
of their absolute values are our estimate of the uncertainties on 
the atmospheric parameters. 

By using this second method we obtained the uncertainties 
AT eff =37 K, Alog(g) =0. 12 and Av, =0.06 km/s, consistent with 
the ones determined with the first method. 

Our best estimate of the total error associated to the abun- 
dance measures is calculated as 



0"tot - -J&EW + <T atm 

listed in the column 7 of Tab. [3] For O cr to t was obtained in a 
different way since its measure was not based on the EW method 
(see Section 6.4). In all the plots for the error bars associated 
with the measure of abundance we adopted <x tot . 

In addition, systematic errors on abundances can be intro- 
duced by systematic errors in the scales of T e ff, log(g), v t , and 
by deviations of the real stellar atmosphere from the model (i.e. 
deviations from LTE). A study of the systematic effects goes be- 
yond the purposes of this study, where we are more interested in 
internal star to star variation in chemical abundances. 

Comparing <x tot with the observed dispersion <x t, s (column 
8) at least for the abundance ratios [Na/Fe], [O/Fe] and [Al/Fe], 
we found <x obs significantly larger than <x tot , suggesting the pres- 
ence of a real spread in the content of these elements within the 
cluster. 



6. The chemical composition of M4 

In the following sections, we will discuss in detail the measured 
chemical abundances. In addition, we will use the abundances of 
C and N from the literature to extend our analysis to those ele- 
ments. The errors we give here are only internal. External ones 
can be estimated by comparison with other works (see Tab.|2]l. 



6. 1 . Iron and iron-peak elements 

The weighted mean of the [Fe/H] found in the 105 cluster mem- 
bers is 

[Fe/H] = -1.07 ± 0.01 dex 

The other Fe-peak elements (Ni and Cr) show roughly the Solar- 
scaled abundance. 

The agreement between cr (, s and cr tot (see Tab. [3]), allows us 
to conclude that there is no evidence for an internal dispersion 
in the iron-peak element content, suggesting that this cluster is 
homogeneous in metallicity down to the <x b s level. 

6.2. a elements 

The chemical abundances for the a elements O, Mg, Si, Ca, and 
Ti are listed in Tab. [2] All the a elements are overabundant. Ca 
and Ti are enhanced by ~0.3 dex, while Si and Mg by ~0.5 dex. 
The results for [O/Fe] will be discussed in detail in Section 6.4. 

Previous investigators (GQO86, BW92 and 199) have also 
found significant overabundances for both Si and Mg. In partic- 
ular, for [Si/Fe] they found values higher than 0.5 dex, slightly 
higher than our results, while our [Mg/Fe] = 0.50 lies in the 
middle of the literature results that range between 0.37 (BW92) 
and 0.68 dex (GQO86). As noted by 199, the abundances of these 
two elements in M4 are higher if compared to the ones found in 
other globular clusters showing a similar metallicity, i. e. M5 
(Ivans et al. 2000). 199 suggest that the higher abundances in 
M4 should be primordial, e.g. due to the chemical composition 
of the primordial site from which the cluster formed. 

For [Ca/Fe] we obtained a scatter of 0.04 dex, differing from 
199, who found <x=0.11 dex (see discussion in Section 6.6). Til 
and Till show the same abundance within 3<x (see Tab. with 
a mean of the two different ionized species of [Ti/Fe] = +0.32 ± 
0.04. 

The average of all the a elements, gives an a enhancement 
for M4 of: 

[or/Fe] = +0.39 + 0.05 dex 

In this case the agreement between cr b s and cr tot again allows us 
to conclude that, down to the <x obs level, there is no evidence of 
internal dispersion for these elements, with the exception of O, 
which we will discuss in Section 6.4. 
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Table 3. Sensitivity of derived UVES abundances to the atmospheric parameters and EWs. We reported the total error due to the 
atmospheric parameters (cr atm ), due to the error in EW measurement (cr EW ), the squared sum of the two (cr tot ), and the observed 
dispersion (<x obs ) for each element. The Fell observed dispersion is not reported because it is necessarily the same as Fel due to the 
method used in this paper to obtain atmospheric parameters. 





AT «(K) 




Avt (km/s) 


u aim 




w tot 


(T u 

XJ obs 




+40 


+0.12 


+0.06 










roi /Fei 












0.04 


0.09 


[Nal/Fe] 


+0.00 


-0.01 


+0.02 


0.02 


0.04 


0.04 


0.17 


[Mgl/Fe] 


-0.01 


-0.01 


+0.01 


0.02 


0.06 


0.06 


0.06 


[All/Fe] 


-0.01 


-0.02 


+0.02 


0.03 


0.06 


0.07 


0.11 


[Sil/Fe] 


-0.04 


+0.03 


+0.02 


0.05 


0.03 


0.06 


0.05 


[Cal/Fe] 


+0.01 


-0.02 


+0.00 


0.02 


0.02 


0.03 


0.04 


[Til/Fe] 


+0.04 


-0.02 


+0.00 


0.04 


0.02 


0.04 


0.05 


[Till/Fe] 


+0.03 


-0.01 


-0.01 


0.03 


0.05 


0.06 


0.06 


[CrI/Fe] 


+0.02 


-0.01 


+0.00 


0.02 


0.06 


0.06 


0.05 


[Fel/H] 


+0.04 


+0.00 


-0.03 


0.05 


0.01 


0.05 


0.05 


[Fell/H] 


-0.03 


+0.06 


-0.01 


0.07 


0.04 


0.08 




[Nil/Fe] 


-0.01 


+0.02 


+0.01 


0.02 


0.02 


0.03 


0.03 


[Ball/Fe] 


+0.05 


-0.04 


-0.03 


0.07 


0.05 


0.09 


0.09 



6.3. Barium 

As in 199 and BW92, we have found a strong overabundance 
of Ba. We have [Ba/Fe] = 0.41 + 0.01 which is about 0.2 dex 
smaller than in 199. In any case, our results confirm that Ba is 
significantly overabundant in M4, at variance with what found 
by GQ086, possibly solving this important issue. As deeply dis- 
cussed by 199, this high Barium content cannot be accounted 
for mixing processes, but must be a signature of the presence 
of s-process elements in the primordial material from which M4 
stars formed. The excess of s-process elements provides some 
support to the idea that the formation of the stars we now ob- 
serve in M4 happened after intermediate-mass AGB stars have 
polluted the environment, or that it lasted for long enough to al- 
low intermediate-mass AGB stars to strongly pollute the lower 
mass forming stars. We do not find any significant dispersion in 
Ba content for the bulk of our target stars, though there are a few 
outliers worth further discussion (see below). 



dances for 24 giant stars from high resolution spectra taken at the 
Lick and McDonald Observatories. With respect to their work, 
our dispersion for [O/Fe] is quite a bit smaller (0.09 dex against 
their value of 0.14 dex), and our average value is larger (see 
Tab. |2). We have seven stars in common with their high reso- 
lution sample, and another seven with their medium resolution 
sample for which they derived only the Oxygen abundance. The 
comparison between their results and the present ones are listed 
in Tab. |H where the last two columns give the differences be- 
tween 199 and our results. We can note that our abundances are 
systematically higher than those of 199. The effect is higher for 
the [O/Fe], but is present for Sodium too, and could reflect the 
difference of ~ 0. 1 dex we found for the iron with respect to 199. 
This systematic difference does not affect the shape of the Na-O 
anticorrelation, but it does emphasize the overabundance of O in 



6.4. Na-O anticorrelation 

Variations in light element abundances is common among GCs, 
and is also present in M4. In particular, Sodium and Oxygen 
have very large dispersions: we obtained o- [Na / Fe ]=0.17 and 
o"[O/Fe]=0.09. A large star to star scatter in Sodium and Oxygen 
abundance in M4 was found also in other previous studies (199, 
D92). 

Sodium and Oxygen show the typical anticorrelation found 
in many other GCs (Gratton et al. 2004), as shown in Fig. [5] 
where the [Na/Fe] values are plotted as a function of [O/Fe]. 
The open blue circles represent the 12 stars for which it was not 
possible to obtain a good estimate of the Oxygen abundance. 

Since the Oxygen abundance was derived from just one spec- 
tral line, we adopted the following the following procedure to de- 
rive a raw estimate of the typical error <x tot associated with each 
measure: we selected the group of stars with [Na/Fe]<0.20, as- 
suming them to be homogeneous in their O content. Under this 
hypothesis, the dispersion in O can be assumed to be due to the 
random error associated with the measured O abundance. We 
obtained o"ro/ Fe j = 0.04 dex, which corresponds to the error bar 
size in Fig.[5]and to the <x tot in Tab. [3] 

The only previous study showing the Na-O anticorrelation in 
M4 was that by 199. They analysed [Na/Fe] and [O/Fe] abun- 
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Fig. 5. [Na/Fe] vs. [O/Fe] abundance ratios. The error bars rep- 
resent the typical errors cr tot from Tab. [3] 
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M4 stars with respect to other clusters with similar metallicity, 
as already noticed by 199: M4 stars formed in an environment 
particularly rich in O and possibly in Na. 

The most interesting result of our investigation is that, thanks 
to the large number of stars in our sample, we can show that the 
[Na/Fe] distribution is bimodal (Fig. [6}. Setting an arbitrary sep- 
aration between the two peaks in Fig. [6] at ~0.2 dex, we obtain, 
for the two groups of stars with higher and lower Na content, a 
mean [Na/Fe] = 0.38 ± 0.01 and [Na/Fe] = 0.07 + 0.01, respec- 
tively. 

In correspondence with the two peaks in the [Na/Fe] 
distribution, there are also two [O/Fe] groups (Fig. [5]): the 
first group is centered at [O/Fe] ~ 0.30, while the second is at 
[O/Fe] ~ 0.47, with a few stars having intermediate Oxygen 
abundance. 

We have followed two methods to trace the distribution of 
the stars on the Na-0 plane. The first method uses the ratio 
[O/Na] between the O and the Na abundances (Carretta et al. 
2006) which appears to be the best indicator to trace the star 
distribution along the Na-O anticorrelation, because this ratio 
continues to vary even at the extreme values along the distri- 
bution. The distribution of the [O/Na] is plotted in Fig. [7] (left 
panel). We can identify two bulks of RGB stars in this plot: 

one at [O/Na] ~ 0.40 and the second one at [O/Na] 0.10. 

There might be another group of stars between [O/Na]~0 and 
[O/Na]~0.25, but we cannot provide conclusive evidence. 

Note that the two groups of stars in the Na-O anticor- 
relation have their corresponding peaks at [O/Fe] ~ 0.47 and 
[O/Fe] ~ 0.30, so both groups are O-enhanced. 
There are no differences in [Fe/H] content between these two 
groups of stars: both of them have [Fe/H] = -1.07 dex. 

We have also analysed the distribution of stars along the Na- 
O anticorrelation with an alternative procedure whose steps are 
briefly described below: 



- A parametric curve following the observed Na-O distribution 
as shown in the inset of the right panel of Fig. [7] has been 
traced; 

- Each observed point in the Na-O distribution has been pro- 
jected on the parametric curve; 

- For each projected point, the distance (P) from the origin of 
the parametric curve (indicated with in the inset of the right 
panel of Fig. [7j has been calculated; 

- A histogram of the P distance distribution has been con- 
structed. 

The histogram is shown in the right panel of Fig. [7] In this case 
also, two peaks are evident: the first one at P ~ 0.28, and the 
second one at P ~ 0.58. 

In conclusion, we obtained the Na-O anticorrelation for a 
large sample (93 RGB stars) in the globular cluster M4. The dis- 
tribution of the objects on the Na-O anticorrelation is clearly bi- 
modal. 

6.5. Mg-AI anticorrelation 

Aluminium abundances have been determined from the EWs of 
the All lines at 6696 A and 6698 A, and Mg abundances from 
the EWs of the Mgl doublet at 6319 A, 6318 A, and of the line 
at 5711 A. Figure [8] shows the [Mg/Fe] as a function of [Al/Fe] 
and [O/Fe] (upper panels) and [Mg/Fe] and [Al/Fe] as a function 
of [Na/Fe] (lower panels). 

In each panel the red line represents the best least squares 
fit to the data. The value a is the slope of the best fit straight 
line y = ax + b. There is no clear Mg-AI anticorrelation (up- 
per left panel of Fig. [8}, although Al is more spread out than 
Mg (o-[ A1/Fe] =0.11 vs. o~ [Mg/Fe] = 0.06). Both Al and Mg are 
overabundant, with average value of [Al/Fe] = 0.54 ± 0.01 and 
[Mg/Fe] = 0.50 + 0.01 (internal errors). 
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Fig. 7. Left: Distribution of stars along the Na-O anticorrelation 
represented by the [O/Na] ratios. Right: Distribution of the pro- 
jected position P of stars on the parametric curve plotted in the 
Fig. 6. Histogram of the distribution of the [Na/Fe] abundances. i nse t pannel. The coordinates of the edges of the curve are indi- 
The error bars represent the Poisson errors. cated. 
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Table 4. Comparison between our [Na/Fe] and [O/Fe] and those of 199 for the stars in common. The Oxygen values were compared 
both with the high resolution sample of 199 ([0/Fe]^ 9 ) and with the medium resolution sample ([O/Fe]^). The last three columns 
give the differences [el/Fe] this work - [el/FeJeg. 



ID 


rOI/Felthk wnrlt 
L i JiniS WOIK 


roi/Fei!;;!, 

L " / Jjyy 


[OI/Fel!" 

L v 1 Jjyy 


INal/Felthk wnrit 

L '"J-/ " w JttlIS WOIK 


rNai/Fei!;;' 

l ' *■"-/ J Jyy 


A(0 (a) ) 


A(0 <b) ) 


A(Na) 


L1411 


+0.28 


+0.20 


+0.07 


+0.51 


+0.43 


+0.08 


+0.21 


+0.08 


L1501 


+0.34 


+0.10 




+0.51 


+0.42 


+0.24 




+0.09 


L1514 


+0.50 


+0.41 


+0.41 


+0.17 


+0.01 


+0.09 


+0.09 


+0.16 


L2519 


+0.48 


+0.37 


+0.37 


+0.02 


-0.19 


+0.11 


+0.11 


+0.21 


L2617 


+0.25 


+0.01 


+0.07 


+0.44 


+0.50 


+0.24 


+0.18 


-0.06 


L3612 


+0.25 


+0.10 


+0.12 


+0.46 


+0.47 


+0.15 


+0.13 


-0.01 


L3624 


+0.48 


+0.29 


+0.27 


+0.18 


+0.10 


+0.19 


+0.21 


+0.08 


L1403 


+0.31 




+0.22 


+0.53 






+0.09 




L2410 


+0.39 




+0.32 


+0.11 






+0.07 




L2608 


+0.43 




+0.37 


+0.24 






+0.06 




L1617 


+0.24 




+0.17 


+0.46 






+0.07 




L4415 


+0.44 




+0.37 


+0.37 






+0.07 




L4421 


+0.28 




+0.07 


+0.44 






+0.21 




L4413 


+0.22 




+0.27 


+0.48 






-0.05 




mean 


+0.35 


+0.21 


+0.24 


+0.35 


+0.25 


+0.16 


+0.11 


+0.08 



The right-bottom panel of Fig. [8] shows that there is a corre- 
lation between [Al/Fe] and [Na/Fe] (higher Al abundances for 
higher Na content), though it is less pronounced than the Na- 
O anticorrelation. Aluminium spans over a range of ~ 0.4 dex, 
while Sodium covers a range of ~ 0.6 dex. The spread in Al con- 
tent of M4 is smaller than in other clusters of similar metallicity, 
like M5 (Ivans et al. 1999) and M13 (Johnson et al. 2005). 
Figure [9] shows that there might also be a correlation between 
the [Al/Fe] and [Ba/Fe]. 

Mg does not clearly correlate neither with Na (left-bottom 
panel of Fig. [8]l nor O (upper-right panel of Fig. [8]). 

The absence of a Mg-Al anticorrelation is somehow surpris- 
ing, in view of the presence of a strong double peaked Na-O 
anticorrelation, and of a correlation between Al and Na and pos- 
sibly Al and Ba. All these correlations seems to indicate the 
presence of material that has gone through the s-process phe- 
nomenon. Assuming that the Na enhancement comes from the 
proton capture process at the expense of Ne, we can also expect 
that a similar process is the basis of the Al enhancement at the 
expense of Mg. Still, we do not find any significant dispersion 
in the Mg content. This problem has already been noted by 199, 
who suggested that the required Mg destruction is too low to be 
measured (see their Section 4.2.2 for further details). 

Here we only add that, despite our large sample, neither 
in the Al nor in the Mg (and Ba) distribution can we see 
any evidence of the two peaks so clearly visible as in the Na 
and O distribution. But, in spite of the lack of a clear Mg- 
Al anticorreletion, we found that Na poor stars have on aver- 
age higher Mg and lower Al. The difference between the me- 
dian Mg and Al abundances for the Na-poor and Na-rich sam- 
ples are: [Mg/Fe] Na _ rich - [Mg/Fe] Na _ poor = -0.04 + 0.01 and 
[Al/Fe] Na _ rich - [Al/Fe] Na _ poor = +0.08 + 0.02 This should be 
consistent with the scenario proposed by 199 (see their Sec. 4.2.2 
for more details) who predicted that a drop of only 0.05 dex in 
Mg is needed to account for the observed increase in the abun- 
dance of Al. (Obviously this requires the further hypothesis that 
all the Na enhanced stars are also Al enhanced.) 



By studying 4 RGB stars from high resolution spectra, D92 
found higher Sodium abundances in the two CN-strong stars 
with respect to the two CN-weak ones (Suntzeff & Smith 1991). 
199 found that the Oxygen abundance is antic orrelated with 
Nitrogen, whereas both Na and Al are more abundant in CN- 
strong giants than in CN-weak ones. 

In our study, two peaks in the Na distribution are evident (see 
Fig. [6}, while there is no equally strong pattern in the Al distri- 
bution (right bottom panel of Fig. [§). The CN band strengths 
of red giants in M4 have been measured by N81, Suntzeff & 
Smith (1991), and 199. Smith & Briley 2005 (SB05) have ho- 
mogenized all the available data for CN band strengths in terms 
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6.6. CN bimodality 

In this section we discuss the correlation of the /13883 CN ab- 
sorption band strength with Na, O, and Al abundances. Previous 
studies revealed a trend of Na and Al with the CN band strength. 



Fig. 8. Bottom panels: [Al/Fe] and [Mg/Fe] abundance ratios 
are plotted as a function of [Na/Fe]. Upper panels: [Mg/Fe] ratio 
as a function of [O/Fe] and [Al/Fe]. In each panel the red line is 
the best least squares fit and a is the slope of this line. The error 
bars represent the typical errors cr tot from Tab. [3] 
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Fig. 9. [Al/Fe] vs. [Ba/Fe] abundance ratios. The error bars rep- 
resent the typical errors cr tot from Tab. [3] 

of the S(3839) index, e. g. the ratio of the flux intensities of the 
Cyanogen band near 3839 A and the nearby continuum. In the 
following, we will use the S(3839) index calculated by SB05. 

In Tab.|5]the CN-index S(3839), the [C/Fe] and the [N/Fe] 
values (taken from SB05 and 199, respectively), and the identi- 
fication as CN-strong (S), -weak (W), and -intermediate (I), are 
listed for our target stars with available CN data (34 stars in to- 
tal). The IDs come from Lee (1977). 

For these targets, Table [6] lists the mean values of the chem- 
ical abundances for CN-S and CN-W(+I) stars. The Na, O, Mg, 
Al and Ca abundances are from this work, the [C/Fe] values 
are given by SB05, while the [N/Fe] are the mean values taken 
from 199 for all their CN-S and CN-W stars, including stars not 
in common with our sample (there are only five stars in our sam- 
ple with nitrogen content available). In the last column the dif- 
ferences A(S - W) between the mean values for CN-S and CN- 
W(+I) stars are listed. 

Fig. [TO] shows the [Na/Fe], [Al/Fe], [O/Fe], [Mg/Fe] and 
[Ca/Fe] values obtained in this work as a function of S(3839). In 
this figure CN-S objects are represented by filled circles, CN-W 
by open circles, and CN-I by open squares. The dotted lines rep- 
resent the mean values of the abundances for the CN-S and CN- 
W(+I) stars; the CN-I stars were included in the CN-W group be- 
cause of their similar behaviour as is clear from Fig. [10] The red 
lines represent the best least squares fit of the data, and a is the 
slope of the line. We note that the CN-S objects clearly show sig- 
nificantly enhanced Na abundances, while the CN-W(+I) stars 
have a lower Na content. The three CN-I stars show low [Na/Fe] 
values typical of CN-weak stars. There is also a systematic dif- 
ference in Al between CN-weak and CN-strong stars, with CN- 
strong stars richer in Al content. 

The CN-W(+I) stars have a higher O content than the CN- 
S as shown in the upper middle panel of Fig. [10] Fig. [TUl shows 
some difference in [Mg/Fe] between the CN-W(+I) and the CN- 
S stars, with CN-W(+I) stars having slightly higher Mg abun- 
dances than the CN-S ones, but this difference is not statistically 
significant (see Table[6]l. 

199 and D92 also found that the scatter in Ca abundance cor- 
relates with the CN strength. 199 found a difference in Ca content 



between CN-strong and CN-weak stars by 0.08 + 0.11 dex. To 
compare, we have also analysed the Ca abundance as a function 
of the CN strength. In the bottom-right panel of Fig.[l0] [Ca/Fe] 
is plotted as a function of S(3839); we can see that the differ- 
ence between the two groups is 0.02 dex, smaller than what was 
found by 199 and not statistically significant (see Table[6]l. 

7. Search for evolutionary effects 

7. 1 . Abundances along the RGB 

According to the results discussed in the previous section, we 
can define two groups of stars: the Na-rich stars, i. e. those 
with [Na/Fe] > 0.2 and low Oxygen content which must be 
associated to the CN-S stars (Fig. [101 . and the Na-poor stars 
with higher Oxygen, i. e. those with [Na/Fe] < 0.2, which cor- 
respond to the CN-W group. It is instructive to look at the po- 
sition of the stars belonging to these two Na groups on the U 
vs. (U - B) CMD and on the color-color (U - B) vs. (B - K) 
diagram (Fig.[TT1and Fig. [12] respectively). The U and B magni- 
tudes come from our WFI photometry, the K magnitude from the 
2MASS catalogue. The Na-rich stars are represented by red tri- 
angles, while the Na-poor stars by blue circles. The CMD has 
been corrected for differential reddening (following the same 
procedure used in Sarajedini et al. 2007) and proper motions. 
It is clear from Fig. [TT] that our sample of stars mostly includes 
RGB stars with only three or four probable AGB stars. 

Interestingly enough, the two groups of Na-rich and Na-poor 
stars form two distinct branches on the RGB. Na-rich stars define 
a narrow sequence on the red side of the RGB, while the Na-poor 
sample populates the blue, more spread out portion of the RGB. 
Even more interestingly, the anomalous broadening of the RGB 
is visible down to the base of the RGB, at U ~ 17.5, indicating 
that the two abundance groups are present all over the RGB, even 
well below the RGB bump, where no deep mixing is expected. 
This evidence further strengthens the idea that the bimodal Na, 
O, CN distribution must have been present in the material from 
which the stars we presently observe in M4 originated. 

FigureQT]unequivocally shows that the Na (CN) dichotomy 
is associated to a dichotomy in the color of the RGB stars. In or- 
der to quantify this split, we calculated for the Na-poor group of 
stars the mean difference in (U - B) color with respect to a fidu- 
cial line representative of the RGB Na-rich stars (see Fig. [T3l l. 
The fiducial line has been obtained as follows: 

- the CMD has been divided into bins of 0.30 mag in U, and a 
median color (U - B) and U magnitude have been computed 
for the Na-rich ([Na/Fe] > 0.2) stars in each bin. The median 
points have been fitted with a spline function which represent 
a first, raw fiducial line; 

- for each Na-rich star, the difference in color with respect to 
the fiducial was calculated, and the 68.27th percentile of the 
absolute values of the color differences, was taken as an es- 
timate of the color dispersion (cr). All stars with a color dis- 
tance from the fiducial larger than 3<x were rejected; 

- the median colors and magnitudes and the cr of each bin were 
redetermined by using the remaining stars. 

The leftmost panel of Fig.[T3"1shows the original CMD, while 
the rightmost one shows the CMD after subtracting from each 
Na-poor star the fiducial line color appropriate for its U magni- 
tude. The differences are indicated as A(U - B). We calculated 
the average A(U - B) of Na-poor stars with U>14.8. This cut 
in magnitude has been imposed by the poor statistics and the 
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Table 5. CN band strengths and Carbon and Nitrogen abundances for our target stars with C, N and CN measurements in the 
literature. We use the star identifications from Lee (1977). In the last column, S refers to CN-strong, W CN-weak and I CN- 
intermediate stars. The C and N abundances are from SB05 and 199, respectively. 



ID S(3839) CN [CI/Fe] [NI/Fe] ID S(3839) CN [CI/Fe] [NI/Fe] 



L1411 


0.54 


s 






1 2626 


0.66 


s 






LI 403 


0.655 


s 


-0.67 




L3705 


0.20 


w 






L1501 


0.61 


s 






L3706 


0.16 


w 






L2422 


0.69 


s 


-0.83 




L2519 


0.39 


I 


-0.50 


+0.22 


L4404 


0.505 


s 


-0.35 




L2623 


0.245 


w 


-0.70 




L1514 




w 


— 1.12 


TU.77 


L3721 


U.JO 


Q 






L4415 


0.53 


s 


-0.67 




L2617 


0.63 


s 


-0.77 


+ 1.05 


L4416 


0.62 


s 


-0.89 




L2711 


0.65 


s 






L4509 


0.62 


s 


-0.44 




L3730 


0.24 


w 






L1512 


0.61 


s 






L3612 


0.65 


s 


-0.73 


+1.02 


L1617 


0.63 


s 


-0.93 




L3621 


0.23 


w 






L1619 


0.62 


s 






L2608 


0.57 


s 


-0.52 




LI 608 


0.60 


s 


-0.63 




L3624 


0.40 


I 


-0.67 


+0.25 


LI 605 


0.74 


s 






L2410 


0.23 


w 


-0.59 




L4630 


0.505 


s 


-0.62 




L4413 


0.655 


s 


-0.32 




L3701 


0.39 


I 






L4421 


0.655 


s 


-0.83 




L3404 


0.12 


w 


-0.40 




L4508 


0.14 


w 








Fig. 10. The abundance ratios [Na/Fe], [O/Fe], [Al/Fe], [Mg/Fe] and [Ca/Fe] are plotted as a function of the CN index S(3839). 
The filled circles represent the CN-S stars, the open circles the CN-W and the open squares the CN-I ones. In each panel, a is the 
slope of the straight line of the best least square fit. The error bars represent the typical errors <x tot from Tab. [3] 

Table 6. Mean abundance ratios and relative rms (cr) for the CN-S and CN-W(+I) groups. In the last column the difference between 
CN-S and CN-W(+I) stars are listed. 



[el/Fe] 


CN-S 


O" 


CN-W(+I) 


tr 


A(S- W(+I)) 


[CI/Fe] 


-0.66 ± 0.05 


0.19 


-0.66 ±0.11 


0.25 


+0.00 


[NI/Fe] 


+ 1.08 + 0.10 


0.22 


+0.42 ±0.19 


0.42 


+0.66 


[OI/Fe] 


+0.33 + 0.01 


0.07 


+0.45 ± 0.01 


0.04 


-0.12 


[Nal/Fe] 


+0.41 ±0.02 


0.07 


+0.11 ±0.02 


0.06 


+0.30 


[All/Fe] 


+0.59 ± 0.01 


0.06 


+0.51 ±0.03 


0.09 


+0.08 


[Mgl/Fe] 


+0.52 ± 0.01 


0.04 


+0.54 ± 0.02 


0.06 


-0.02 


[Cal/Fe] 


+0.26 ± 0.01 


0.04 


+0.24 ± 0.01 


0.03 


+0.02 
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Fig. 11. U vs. (U - B) CMD from WFI photometry. The distribution of the Na content for the UVES stars is shown in the inset. 
The stars belonging to the two different Na groups are represented in two different colors: the red triangles represent the stars with 
[Na/Fe] > 0.2 dex, and the blue circles the stars with [Na/Fe] < 0.2 dex. 



presence of the probable AGB stars at brighter magnitudes. The 
mean A(U - B) value is -0. 17 ± 0.02 (dotted line in the left most 
panel of Fig. H~3b. 

The significant difference between the mean colors of the 
2 groups of stars is a further evidence of the presence of two 
different stellar populations in the RGB of M4, to be associated 
with the different content of Na, and, because of the discussed 
correlations to different O, N, and C content. 

In order to understand the origin of the photometric di- 
chotomy, using SPECTRUM, we simulated two synthetic spec- 
tra, one representative of the Na-rich, CN-strong stars and one 
for the Na-poor, CN-weak ones. The two synthetic spectra were 
computed using as atmospheric parameters the mean values 
measured for the sample of our stars for which there are liter- 
ature data on the CN band strengths (see Tab. [5), and assuming 
for the C, N, O, Na abundances the average values calculated 
for the CN-W(+I) and the CN-S groups and listed in Tab.|6] We 
then multiplied the two spectra by the efficency curve of our U 
and B photometric bands (Fig.[l4l lower panels), and, finally, we 
calculated the difference between the resulting fluxes. The upper 
panels of Fig.fPflshow the differences between the two simulated 
spectra as a function of the wavelength for the U (left panel) and 
B (right panel) panel. It is clear that the strength of the CN and 
NH bands strongly influences the U-B color. The NH band 
around 3360 A, and the CN bands around 3590 A, 3883 A and 
4215 A are the main contributors to the effect. 

The differences in magnitude between the CN-S 
and the CN-W(+I) simulated spectra in the two bands 
are: A(U)st r0 ng-Weak = +0.06 and A(B)s tl -ong-Weak = +0.02. 
Consequently, the expected color difference between the two 




B-K 

Fig. 12. (U - B) vs. (B - K) diagram: U and B come from WFI 
photometry, K from 2MASS catalogue. The UVES stars belong- 
ing to the two different Na groups are represented as in Fig.Qj] 



groups of stars is A(U - B)strong-Weak = +0.04. This value goes 
in the same direction, but it is smaller than the observed one 
(A(U - B)st ro ng-Weak = +0.17 + 0.02). However, we note that 
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Fig. 13. (Left panel) Distribution along the RGB of the stars with 
measured abundances. The Na-rich (red) and the Na-poor (blue) 
stars are represented as in Fig. QT| and Fig. Q~2] In the middle 
panel the N re d/Nbi ue ratios for the 4 selected magnitude bins are 
shown with their relative Poisson errors. The rightmost panel 
shows the color difference between each analyzed star and a ref- 
erence fiducial line (dotted line). 



our procedure uses simulated spectra with average abundances, 
and therefore should be considered a rough simulation. We 
cannot exclude the possibility that other effects (perhaps related 
to the structure and evolution of the stars with different chemical 
content, or effects on the stellar atmosphere associated with the 
complex distribution of the chemical abundances in addition 
to the CN band), might further contribute to the photometric 
dichotomy in the RGB. Surely, our simulations show that the 
CN-bimodality affects the (U - B) color and can be at least 
partly responsible for the observed spread in the U vs. (U - B) 
CMD. 

Finally, we investigated whether the chemical and photo- 
metric dichotomy is related to the evolutionary status along the 
RGB. To this end, we have calculated the fractions of stars in the 
2 Na-groups at different magnitudes along the RGB. As shown 
in the middle panel of Fig. Qj] we divided the RGB in 4 mag- 
nitude bins containing the same number of stars with measured 
metal content, and calculated the ratio between the number of 
Na-rich (N re( j) and the number of Na-poor (Nbiue) stars in each 
bin. The N re d/Nbi U e ratios with their associated Poisson errors 
as a function of the magnitude are plotted the middle panel of 
Fig. [O] The N re d/Nbi ue values in the different magnitude bins 
are the same within the errors. A similar result is obtained if 
we divide the RGB in two bins only, but containing a number 
of stars twice as large as in the previous experiment. In this 
case, we have N red /N b i ue = 1.52 ± 0.44 for 14.0 < U < 14.9, 
N re d/Nbiue = 1.82 ± 0.54 for 14.9 < U < 15.8. 

Figure [15] shows the dependence of Na, Al, and S(3839) in- 
dex as a function of the magnitude U (upper panel) and B (lower 
panel). Again, we do not see any trend of the abundances with 
the position of the stars along the RGB. 



In conclusion, there is no evidence for a dependence of the 
Na content and (because of the discussed correlations) of the O, 
Al abundances or of the CN strength with the evolutionary status 
of the stars. 

7.2. The RGB progeny on the HB 

In the previous sections, we have identified two groups of stars 
with distinct Na, O, CN content, and which populate two distinct 
regions in the U vs. U - B CMD. 

We note that Yong & Grundhal (2008), by studying 8 bright 
red giants in the GC NGC 1851 found a large star to star abun- 
dance variation in the Na, O, and Al content, and suggested that 
these abundance anomalies could be associated with the pres- 
ence of the two stellar populations identified by Milone et al. 
(2008), a situation apparently very similar to what we have found 
in M4. Cassisi et al. (2008) suggest that extreme CNO variations 
could be the basis of the SGB split found in NGC 1851. As dis- 
cussed in Milone et al. (2008) there is some evidence of a split of 
the RGB in NGC 1851 similar to what we have found in this pa- 
per for M4. We have investigated all available HST data of M4, 
including CMDs from very high precision photometry based on 
ACS/HST images; we could find no evidence for a SGB split (as 
in NGC 1851) or evidence for a MS split (as in NGC 2808 or 
at Cen). 

It is also useful to investigate where the progeny of the two 
RGB populations is along the HB. Milone et al. (2008) suggest 
that the bimodal HB of NGC 1851 can be interpreted in terms 
of the presence of two distinct stellar populations in this cluster, 
and tentatively link the bright SGB to the red HB (RHB) and the 
fainter SGB to the blue HB (BHB). M4 also has a bimodal HB, 
well populated on both the red and blue side of the RR Lyrae 
gap. Could the HB morphology of M4 be related to the CN bi- 
modality? N81, studying the CN band strengths on a sample of 
45 giant stars in M4, suggested the possibility of a relation be- 
tween the HB morphology and the CN-bimodality. 

To test the possible relation between the chemical and con- 
sequent RGB bimodality of M4 and the morphology of its HB, 
we used the WFI@2.2m instrumental photometry by Anderson 
et al. 2006 (Fig. |T6t corrected for differential reddening, and se- 
lected the cluster members using the measured proper motions. 
We remark here the fact that this photometry is not the one we 
use in the whole paper. The homogenization between the two 
photometries is a difficult and long process (mainly because the 
fact that the B filters for the two datasets are different) beyond 
the scope of this work. 

In our test, we used the V vs. (B - V) CMD, as in this pho- 
tometric plane the two components of the HB are more clearly 
distinguishable. 

In Fig. [16] different symbols show the red and blue HB stars 
that we selected. The ratio between the number (Nbhb) of the 
HB stars bluer than the RR-Lyrae instability strip (BHB), and 
the total HB stars (N TO t) is: 

Nbhb /N T ot = 0.56 + 0.10 

while the ratio between the Na rich stars (Nn 3 r) and the total 
number of stars in our spectroscopic sample (Ns) is: 

N N aR/N s = 0.64 ±0.10 

where the associated uncertainties are the Poisson errors. 
Within the statistical uncertainties, accounting for the different 
evolutionary times along the HB, we can tentatively associate 
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Fig. 15. Dependence of Na, Al, and S(3839) index as a function of the magnitude U (upper panels) and B (lower panels). 



the Na-rich stars with the BHB, and the Na-poor stars with the 
RHB. A direct measurement of the metal content of the HB stars 
is needed in order to support this suggestion.. 

8. Conclusions 

We have presented high resolution spectroscopic analysis of 105 
RGB stars in the GC M4 from UVES data. 

We have found that M4 has an iron content 
[Fe/H] = -1.07 ± 0.01 (the associated error here is the 
internal error only), and an a element overabundance 
[a/Fe] = +0.39 + 0.05. Si and Mg are more overabundant 
than the other a elements, suggesting a primordial overabun- 
dance for these elements. Moreover, we also find a slight 
overabundance of Al. The [Na/Fe] versus [O/Fe] ratios follow 
the well-known Na-O anticorrelation, signature of proton 
capture reactions at high temperature. No Mg-Al anticorrelation 
was found. 



We find a strong dichotomy in Na abundance, and show that 
it must be associated to a CN bimodality. Tab. [7] synthesizes the 
spectroscopic and photometric differences we found for the two 
groups of stars. Comparing our results to those in the literature 
on the CN band strength, the CN-strong stars appear to have 
higher content in Na and Al, and lower O than the CN-weak 
ones. Apparently, M4 hosts two different stellar populations. 
This fact is evident from a chemical abundance distribution, and 
is also confirmed by photometry. In fact, an inspection of the U 
vs. (U - B) CMD reveals a broadened RGB, and, as shown in 
our investigation, the two Na groups of stars occupy different 
positions (have different colors) along the RGB. Our sample of 
stars is composed of objects with U brighter than ~16, but, as 
shown in Fig. QT| the RGB appears to be similarly spread from 
the level of the SGB to the RGB tip. Since our photometry has 
been corrected for differential reddening, we conclude that the 
observed RGB spread at fainter magnitudes must also be due to 
the metal content dichotomy. 
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Table 7. Mean abundance ratios and relative rms (<x) for the Na-rich and Na-poor groups. In the last column, the difference between 
Na-rich and Na-poor stars are listed. 



[el/Fe] 


Na-rich 




Na- 


joor 


a 


A(Na rich 


-poor) 


[OI/Fe] 


+0.34 


±0.01 


0.08 


+0.47 


±0.01 


0.04 


-0.13 ± 


0.01 


[Nal/Fe] 


+0.38 


±0.01 


0.08 


+0.07 


±0.01 


0.06 


+0.31 ± 


0.01 


[Mgl/Fe] 


+0.48 


±0.01 


0.06 


+0.52 


±0.01 


0.05 


-0.04 ± 


0.01 


[All/Fe] 


+0.57 


±0.01 


0.08 


+0.49 


±0.02 


0.13 


+0.08 ± 


0.02 


[Sil/Fe] 


+0.49 


±0.01 


0.04 


+0.46 


±0.01 


0.06 


+0.03 ± 


0.01 


[Cal/Fe] 


+0.28 


±0.01 


0.04 


+0.28 


±0.01 


0.04 


+0.00 ± 


0.01 


[Til/Fe] 


+0.29 


±0.01 


0.03 


+0.29 


±0.01 


0.05 


+0.00 ± 


0.01 


[Till/Fe] 


+0.35 


±0.01 


0.07 


+0.34 


±0.01 


0.05 


+0.01 + 


0.01 


[CrI/Fe] 


-0.04 


±0.01 


0.05 


-0.03 


±0.01 


0.05 


-0.01 + 


0.01 


[Fel/H] 


-1.07 


±0.01 


0.05 


-1.07 


±0.01 


0.05 


+0.00 ± 


0.01 


[Nil/Fe] 


+0.02 


±0.01 


0.03 


+0.01 


±0.01 


0.03 


+0.01 ± 


0.01 


[Ball/Fe] 


+0.42 


±0.01 


0.07 


+0.40 


±0.02 


0.13 


+0.02 ± 


0.02 



fjy [WFI pixel] 
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Fig. 16. In the top panels, the proper motions diagrams show the 
separation of the probable cluster stars (left panel) from the field 
(right panel). In the bottom panels, the cluster and field CMD (V 
vs. B-V) are shown for the stars within the red circle in their 
respective upper panels. In the M4 CMD, the stars on the red and 
on the blue side of the RR Lyrae instability strip are represented 
with red and blue symbols respectively. For this figure, the pho- 
tometry, which is not photometrically calibrated, is taken from 
Anderson et al. (2006). 



We did not find any evident dependence of the chemical 
abundance distribution on the evolutionary status (along the 
RGB) of the target stars, from below the RGB bump to the RGB 
tip. This is an additional indication that the abundance spread 
must be primordial. 

The two groups of stars we have identified both spectroscop- 
ically and photometrically seem to be due to the presence of two 
distinct populations of stars. The abundance anomalies are very 
likely due to primordial variations in the chemical content of the 
material from which M4 stars formed, and not to different evo- 
lutionary paths of the present stellar population of M4. This is 
somehow surprising, because of the relatively small mass of M4 



- it is an order of magnitude smaller than the mass of to Cen, 
NGC 2808, NGC 1851, NGC 6388, the other clusters in which 
a multiple stellar generation has been confirmed. Where did the 
gas which polluted the material for the CN-Na rich stars come 
from? Has it been ejected from a first generation of stars? How 
could it stay within the shallow gravitational potential of M4? 

The idea that M4 hosts two generations of stars makes the 
multipopulation phenomenon in GCs even more puzzling than 
originally thought. It becomes harder and harder to accept the 
idea that the phenomenon can be totally internal to the cluster, 
unless this object is what remains of a much larger system (a 
larger GC or the nucleus of a dwarf galaxy?). Surely, Because 
its orbit involves frequent passages at high inclination through 
the Galactic disk, always at distances from the Galactic center 
smaller than 5 kpc (see Dinescu et al. 1999), M4 must have been 
strongly affected by tidal shocks, and therefore it might have 
been much more massive in the past. 
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Table 8. Coordinates, atmospheric parameters, U,B,V,I C , and 2MASS J,H,K magnitudes for the analised stars. All magnitudes are not corrected 
for differential reddening. 



ID 


RA( U ) 


DEC( U ) 


T et f (K) 


log(g) 


v t (km/; 


19925 


245.81007350 


-26.60153694 


4050 


1.20 


1.67 


20766 


245.80811060 


-26.55680649 


4400 


1.80 


1.45 


21191 


245.81900040 


-26.53602484 


4270 


1.60 


1.60 


21728 


245.81082240 


-26.51303931 


4525 


2.00 


1.42 


22089 


245.82041320 


-26.49663964 


4700 


2.28 


1.36 


24590 


245.91667710 


-26.65088182 


4850 


2.66 


1.35 


25709 


245.92048770 


-26.62319472 


4680 


2.20 


1.38 


26471 


245.88130650 


-26.60513670 


4800 


2.40 


1.28 


26794 


245.90940720 


-26.59866016 


4800 


2.45 


1.44 


27448 


245.95330690 


-26.58586889 


4310 


1.57 


1.58 


28103 


245.84921830 


-26.57493123 


3860 


0.50 


1.62 


28356 


245.97324550 


-26.57086633 


4600 


2.22 


1.53 


28707 


245.93837450 


-26.56579907 


4880 


2.74 


1.31 


28797 


245.97434580 


-26.56458317 


4640 


2.35 


1.36 


28847 


245.90471930 


-26.56390736 


4780 


2.40 


1.27 


28977 


245.94551720 


-26.56218478 


4680 


2.33 


1.40 


29027 


245.84810750 


-26.56146317 


4720 


2.40 


1.41 


29065 


245.94135950 


-26.56081120 


4650 


2.10 


1.41 


29171 


245.89969150 


-26.55958664 


4880 


2.64 


1.26 


29222 


245.84606690 


-26.55896241 


4720 


2.50 


1.36 


29272 


245.86286750 


-26.55825152 


4780 


2.50 


1.26 


29282 


245.85544600 


-26.55811142 


4650 


2.30 


1.42 


29397 


245.87588880 


-26.55671383 


4600 


1.50 


1.78 


29545 


245.92749790 


-26.55485385 


4880 


2.61 


1.18 


29598 


245.88908940 


-26.55414514 


4840 


2.50 


1.40 


29693 


245.88094110 


-26.55307864 


4360 


1.10 


1.88 


29848 


245.89219050 


-26.55108648 


4780 


2.52 


1.24 


30209 


245.89243520 


-26.54681482 


4880 


2.62 


1.32 


30345 


245.94994400 


-26.54524113 


4850 


2.73 


1.31 


30450 


245.86144840 


-26.54411333 


4760 


2.53 


1.35 


30452 


245.88768830 


-26.54411091 


4830 


2.56 


1.25 


30549 


245.88815200 


-26.54294543 


4830 


2.52 


1.28 


30598 


245.89661740 


-26.54235225 


4360 


1.75 


1.47 


30653 


245.88793300 


-26.54184062 


4660 


2.30 


1.25 


30675 


245.90236690 


-26.54157800 


4830 


2.58 


1.35 


30711 


245.85156170 


-26.54123206 


4560 


2.25 


1.46 


30719 


245.90310860 


-26.54116153 


4810 


2.65 


1.24 


30751 


245.94302430 


-26.54089574 


4430 


1.78 


1.47 


30924 


245.89624320 


-26.53883396 


4810 


2.60 


1.28 


30933 


245.95895830 


-26.53873652 


4800 


2.63 


1.30 


31015 


245.89710880 


-26.53795346 


4800 


2.47 


1.37 


31306 


245.92607560 


-26.53505913 


4900 


2.87 


1.33 


31376 


245.91219110 


-26.53432365 


4800 


2.59 


1.36 


31532 


245.89561610 


-26.53290071 


4770 


2.60 


1.21 


31665 


245.93007130 


-26.53164468 


4650 


2.17 


1.34 


31803 


245.90925770 


-26.53023616 


4850 


2.60 


1.34 


31845 


245.85370350 


-26.52982260 


4700 


2.42 


1.31 


32055 


245.91709840 


-26.52773869 


4300 


1.57 


1.52 


32121 


245.83812920 


-26.52713080 


4840 


2.58 


1.33 


32151 


245.83708630 


-26.52688433 


4770 


2.43 


1.38 


32317 


245.86225000 


-26.52536111 


4510 


1.88 


1.43 


32347 


245.92079170 


-26.52505556 


4640 


2.22 


1.38 


32583 


245.89985980 


-26.52267062 


4850 


2.54 


1.34 


32627 


245.97840770 


-26.52219757 


4750 


2.42 


1.30 


32700 


245.91714040 


-26.52145861 


4560 


2.12 


1.36 


32724 


245.92395710 


-26.52127838 


4850 


2.73 


1.30 


32782 


245.89106390 


-26.52065353 


4880 


2.60 


1.24 


32871 


245.85272450 


-26.51984366 


4770 


2.48 


1.24 


32874 


245.91452240 


-26.51985129 


4600 


2.04 


1.38 


32933 


245.83915100 


-26.51926550 


4430 


1.42 


1.78 



u 



B 



I f 



H 



K 



14.83 
14.89 
14.79 
15.05 
15.28 
16.89 
15.11 
15.38 
15.69 

14.75 
14.93 
15.59 
15.13 
15.43 
15.22 
15.18 
15.08 
15.72 
15.28 
15.32 
15.17 
14.60 
15.53 
15.65 
14.67 
15.48 
15.52 
15.64 
15.44 
15.72 
15.52 
14.72 
15.23 
15.52 
15.04 
15.79 
14.83 
15.48 
15.46 
14.91 
15.79 
15.37 
15.28 
14.94 
15.77 
15.23 
14.71 
15.49 
15.41 
15.06 



12.76 
13.53 
13.21 
13.86 
14.36 
14.67 
14.24 
14.70 
14.94 
13.24 
12.52 
13.92 
14.84 
14.22 
14.76 
14.31 
14.40 
14.27 
14.95 
14.42 
14.64 
14.17 
13.53 
14.96 
14.89 
13.27 
14.75 
14.87 
15.00 
14.59 
15.07 
14.86 
13.49 
14.40 
14.77 
13.86 
15.13 
13.61 
14.79 
14.70 
14.21 
15.12 
14.52 
14.53 
14.12 
15.14 
14.34 
13.36 
14.67 
14.59 
13.93 



11.04 
12.10 
11.70 
12.51 
13.14 
13.44 
12.96 
13.49 
13.71 
11.73 
10.71 
12.57 
13.61 
12.93 
13.54 
13.02 
13.16 
12.98 
13.73 
13.17 
13.42 
12.88 
12.20 
13.79 
13.68 
11.80 
13.52 
13.66 
13.84 
13.34 
13.88 
13.66 
12.05 
13.14 
13.58 
12.49 
13.97 
12.20 
13.61 
13.49 
13.00 
13.97 
13.29 
13.33 
12.85 
13.99 
13.09 
11.88 
13.45 
13.38 
12.56 



8.97 

10.33 

9.85 

10.78 

11.54 

11.88 

11.31 

11.87 

12.08 

9.86 

8.51 

10.87 

11.99 

11.30 

11.88 

11.36 

11.48 

11.31 

12.09 

11.52 

11.76 

11.18 

10.47 

12.23 

12.05 

9.96 

11.87 

12.05 

12.31 

11.69 

12.26 

12.05 

10.23 

11.50 

12.00 

10.73 

12.43 

10.45 

12.03 

11.92 

11.38 

12.44 

11.68 

11.76 

11.18 

12.46 

11.44 

10.06 

11.84 

11.75 



15.46 


14.80 


13.61 


12.03 


15.01 


14.34 


13.14 


11.59 


14.88 


13.86 


12.52 


10.80 


15.53 


14.82 


13.62 


12.04 


15.81 


15.10 


13.94 


12.38 


15.65 


14.92 


13.72 


12.09 


14.82 


13.99 


12.72 


11.07 


14.54 


13.29 


11.87 


10.08 



7.59 

9.07 

8.50 

9.50 

10.34 

10.79 

10.16 

10.73 

10.92 

8.57 

7.01 

9.67 

10.80 

10.12 

10.67 

10.14 

10.29 

10.10 

10.88 

10.32 

10.54 

9.91 

9.25 

11.08 

10.86 

8.63 

10.69 

10.84 

11.19 

10.47 

11.07 

10.76 

8.83 

10.29 



6.63 

8.24 

7.67 

8.75 

9.61 

10.11 

9.45 

10.05 

10.25 

7.71 

6.03 

8.92 

10.14 

9.43 

9.98 

9.44 

9.58 

9.42 

10.16 

9.59 

9.86 

9.20 

8.54 

10.43 

10.19 

7.79 

10.00 

10.15 

10.55 

9.78 

10.38 

10.13 

8.03 

9.56 



6.38 

8.07 

7.45 

8.51 

9.47 

9.94 

9.28 

9.90 

10.07 

7.53 

5.75 

8.72 

9.95 

9.21 

9.77 

9.26 

9.39 

9.24 

10.00 

9.41 

9.67 

8.98 

8.32 

10.24 

10.01 

7.64 

9.83 

9.99 

10.39 

9.59 

10.22 

9.99 

7.75 

9.38 



9.44 8.66 8.46 



9.19 


8.36 


8.16 


9.55 


8.81 


8.40 


10.77 


10.12 


9.91 


10.17 


9.48 


9.29 


11.32 


10.64 


10.50 


10.46 


9.79 


9.59 


10.54 


9.85 


9.72 


9.96 


9.25 


9.04 


11.31 


10.64 


10.51 


10.16 


9.49 


9.29 


8.71 


7.90 


7.66 


10.64 


9.93 


9.75 


10.54 


9.84 


9.70 


9.51 


8.76 


8.57 


9.92 


9.20 


9.00 


10.80 


10.16 


9.95 


10.48 


9.79 


9.61 


9.54 


8.81 


8.57 


10.88 


10.20 


10.04 


10.88 


10.22 


10.01 


9.86 


9.13 


8.93 


8.81 


8.02 


7.79 
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ID 


RA( U ) 


DEC( U ) 


T eff (K) 


log(g) 


v, (km/s) 


U 


B 


V 


Ic 


J 


H 


K 


32968 


245.89176500 


-26.51884367 


4630 


2.17 


1.30 


15.06 


14.21 


12.96 


11.32 


10.07 


9.41 


9.15 


32988 


245.90955560 


-26.51865553 


4850 


2.63 


1.22 


15.43 


14.88 


13.75 


12.22 


11.08 


10.40 


10.22 


33069 


245.94851290 


-26.51783666 


4940 


3.05 


1.36 


15.45 


14.71 


13.53 


11.99 


10.86 


10.19 


10.02 


33195 


245.86426280 


-26.51654913 


4620 


2.38 


1.26 


15.30 


14.39 


13.12 


11.45 


10.23 


9.50 


9.36 


33414 


245.84207930 


-26.51444711 


4840 


2.51 


1.28 


15.51 


14.84 


13.67 


12.08 


10.90 


10.21 


10.09 


33617 


245.86252960 


-26.51233938 


4720 


2.35 


1.26 


15.38 


14.57 


13.36 


11.76 


10.59 


9.86 


9.68 


33629 


245.83816250 


-26.51221901 


4930 


2.80 


1.33 


15.77 


15.06 


13.90 


12.34 


11.17 


10.48 


10.36 


33683 


245.89914420 


-26.51176320 


4800 


2.57 


1.18 


15.29 


14.66 


13.48 


11.91 


10.78 


10.09 


9.89 


33788 


245.91070220 


-26.51060047 


4700 


2.37 


1.33 


15.15 


14.24 


12.99 


11.38 


10.16 


9.44 


9.27 


33900 


245.84591550 


-26.50933155 


4770 


2.48 


1.27 


15.48 


14.69 


13.51 


11.93 


10.67 


10.03 


9.92 


33946 


245.84147130 


-26.50887807 


4800 


2.62 


1.15 


15.63 


14.89 


13.71 


12.14 


10.95 


10.29 


10.15 


34006 


245.87135820 


-26.50826984 


4320 


1.67 


1.61 


14.92 


13.37 


11.87 


9.97 


8.59 


7.81 


7.51 


34130 


245.88104790 


-26.50687941 


4550 


2.08 


1.40 


14.98 


13.87 


12.54 


10.82 


9.58 


8.80 


8.60 


34167 


245.90419020 


-26.50638230 


4950 


2.60 


1.40 


15.63 


15.08 


13.95 


12.44 


11.31 


10.67 


10.52 


34240 


245.90241860 


-26.50557893 


4470 


1.95 


1.41 


14.86 


13.74 


12.37 


10.63 


9.35 


8.62 


8.43 


34502 


245.92046770 


-26.50260567 


4860 


2.70 


1.33 


15.49 


14.78 


13.62 


12.09 


10.96 


10.26 


10.15 


34579 


245.93359340 


-26.50170174 


4330 


1.59 


1.49 


14.67 


13.27 


11.83 


10.07 


8.77 


7.91 


7.71 


34726 


245.85009660 


-26.50013473 


4600 


2.24 


1.35 


15.17 


14.14 


12.86 


11.16 


9.89 


9.13 


8.99 


35022 


245.94949140 


-26.49673818 


4850 


2.51 


1.36 


15.50 


14.79 


13.62 


12.10 


10.96 


10.31 


10.11 


35061 


245.92398380 


-26.49632071 


4860 


2.67 


1.23 


15.68 


15.02 


13.85 


12.31 


11.15 


10.48 


10.34 


35455 


245.92832620 


-26.49099726 


4600 


2.10 


1.29 


14.89 


13.99 


12.71 


11.05 


9.83 


9.08 


8.92 


35487 


245.95940170 


-26.49051723 


4850 


2.67 


1.24 


15.33 


14.61 


13.45 


11.95 


10.82 


10.18 


9.99 


35508 


245.85414850 


-26.49026951 


4780 


2.48 


1.18 


15.66 


14.96 


13.80 


12.23 


11.05 


10.37 


10.18 


35571 


245.88484770 


-26.48951670 


4880 


2.79 


1.10 


15.66 


15.07 


13.93 


12.39 


11.25 


10.61 


10.42 


35627 


245.95411440 


-26.48867681 


4830 


2.37 


1.20 


15.48 


14.94 


13.80 


12.30 


11.14 


10.50 


10.34 


35688 


245.91086870 


-26.48777246 


4720 


2.25 


1.33 


15.32 


14.48 


13.27 


11.68 


10.46 


9.75 


9.57 


35774 


245.87296550 


-26.48657771 


4450 


1.92 


1.44 


14.88 


13.68 


12.32 


10.59 


9.32 


8.57 


8.35 


36215 


245.91658540 


-26.48028356 


4300 


1.59 


1.53 


14.74 


13.29 


11.80 


9.96 


8.61 


7.82 


7.57 


36356 


245.90812220 


-26.47815339 


4820 


2.66 


1.26 


15.52 


14.76 


13.57 


12.01 


10.82 


10.17 


9.97 


36929 


245.92953710 


-26.46873630 


4820 


2.55 


1.28 


15.50 


14.75 


13.59 


12.07 


10.93 


10.25 


10.07 


36942 


245.97083610 


-26.46848900 


4800 


2.66 


1.23 


15.43 


14.74 


13.56 


12.04 


10.90 


10.23 


10.08 


37215 


245.89085360 


-26.46381984 


4790 


2.50 


1.21 


15.61 


14.90 


13.72 


12.15 


10.94 


10.26 


10.11 


38075 


245.84038430 


-26.44634617 


4800 


2.54 


1.25 


15.61 


14.83 


13.65 


12.05 


10.82 


10.18 


9.99 


38383 


245.84437580 


-26.43947709 


4590 


1.87 


1.45 


15.02 


14.03 


12.74 


11.02 


9.74 


8.99 


8.79 


38399 


245.88487200 


-26.43904923 


4730 


2.35 


1.38 


15.23 


14.37 


13.14 


11.51 


10.29 


9.56 


9.39 


38896 


245.91189860 


-26.42854610 


4760 


2.53 


1.31 


15.26 


14.36 


13.14 


11.55 


10.33 


9.64 


9.46 


42490 


246.00378460 


-26.59808973 


4570 


2.08 


1.41 


14.95 


13.94 


12.56 


10.84 


9.71 


8.95 


8.76 


42620 


245.99313500 


-26.59417329 


4600 


2.05 


1.37 


14.84 


13.91 


12.59 


10.94 


9.85 


9.11 


8.97 


43370 


246.03196790 


-26.57365632 


4920 


2.80 


1.32 


15.63 


15.03 


13.85 


12.31 


11.29 


10.66 


10.50 


44243 


246.02382800 


-26.55029031 


4860 


2.80 


1.17 


15.51 


14.91 


13.75 


12.24 


11.23 


10.59 


10.42 


44595 


246.06621120 


-26.54096589 


4750 


2.40 


1.40 


15.02 


14.23 


13.00 


11.45 


10.42 


9.75 


9.58 


44616 


245.98944100 


-26.54022147 


4620 


2.20 


1.44 


14.99 


13.99 


12.65 


10.99 


9.82 


9.11 


8.89 


45163 


246.00502890 


-26.52572879 


4770 


2.40 


1.26 


15.27 


14.59 


13.40 


11.88 


10.81 


10.17 


9.98 


45895 


246.01787760 


-26.50775258 


4720 


2.25 


1.34 


14.93 


14.22 


12.99 


11.43 


10.38 


9.71 


9.52 


5359 


245.91226500 


-26.36929347 


4800 


2.44 


1.28 


15.60 


14.78 


13.58 


12.01 


10.87 


10.21 


10.05 
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ID 


[Fel/H] 


[OI/Fe] 


[Nal/FeW 


[Mgl/Felmj-E 


[All/Fe] 


[Sil/Fe] 


[Cal/Fe] 


[Til/Fe] 


[Tffl/Fe] 


[CrI/Fe] 


[Nil/Fe] 


[BaII/Fe] 


19925 


-1.02 


0.28 


0.51 


0.43 


0.68 


0.50 


0.19 


0.48 


0.29 


-0.04 


0.09 


0.64 


20766 


-1.05 


0.31 


0.53 


0.49 


0.62 


0.49 


0.26 


0.31 


0.39 


-0.01 


0.04 


0.45 


21191 


-1.06 


0.34 


0.51 


0.55 


0.59 


0.54 


0.24 


0.31 


0.32 


-0.00 


0.03 


0.53 


21728 


-1.06 


0.31 


0.37 


0.52 


0.68 


0.47 


0.28 


0.32 


0.30 


-0.04 


0.03 


0.48 


22089 


-1.06 


0.35 


0.50 


0.30 


— 


0.53 


0.31 


0.25 


0.32 


-0.15 


0.03 


0.46 


24590 


-1.07 


0.50 


0.30 


0.61 


— 


0.50 


0.27 


0.29 


0.31 


-0.09 


0.04 


0.33 


25709 


-1.13 


0.44 


0.34 


0.54 


0.59 


0.53 


0.35 


0.26 


0.33 


-0.02 


-0.01 


0.35 


26471 


-1.10 


0.41 


0.09 


0.59 


0.58 


0.45 


0.30 


0.23 


0.31 


0.00 


-0.01 


0.28 


26794 


-1.17 


— 


0.36 


0.53 


— 


0.50 


0.34 


0.34 


0.39 


-0.01 


0.03 


0.33 


27448 


-1.12 


0.51 


0.11 


0.50 


0.42 


0.45 


0.29 


0.39 


0.36 


-0.02 


0.03 


0.35 


28103 


-1.08 


0.50 


0.17 


0.41 


0.33 


0.55 


0.16 


0.48 


0.26 


-0.11 


0.01 


0.26 


28356 


-1.14 


0.44 


0.37 


0.55 


0.57 


0.53 


0.28 


0.29 


0.37 


0.04 


0.05 


0.39 


28707 


-1.03 


— 


0.22 


0.40 


— 


0.43 


0.20 


0.27 


0.26 


-0.06 


0.02 


0.45 


28797 


-1.12 


0.35 


0.44 


0.55 


0.57 


0.57 


0.26 


0.28 


0.44 


-0.12 


0.01 


0.53 


28847 


-1.16 


0.52 


0.08 


0.52 


— 


0.51 


0.32 


0.31 


0.30 


-0.09 


-0.04 


0.46 


28977 


-1.14 


0.39 


0.40 


0.51 


0.65 


0.57 


0.30 


0.31 


0.35 


-0.02 


0.05 


0.42 


29027 


-1.10 


0.51 


0.02 


0.48 


— 


0.44 


0.27 


0.27 


0.35 


-0.01 


-0.03 


0.34 


29065 


-1.12 


0.45 


0.17 


0.61 


— 


0.48 


0.27 


0.21 


0.35 


-0.06 


-0.01 


0.35 


29171 


-0.99 


0.22 


0.41 


0.43 


— 


0.49 


0.32 


0.32 


0.31 


-0.04 


-0.04 


0.36 


29222 


-1.04 


0.36 


0.24 


0.44 


0.44 


0.47 


0.25 


0.27 


0.58 


-0.06 


0.05 


0.39 


29272 


-1.11 


0.53 


0.05 


0.47 


— 


0.43 


0.25 


0.30 


0.33 


0.01 


0.02 


0.35 


29282 


-1.06 


0.28 


0.42 


0.50 


0.61 


0.49 


0.25 


0.23 


0.34 


-0.07 


0.02 


0.27 


29397 


-1.12 


0.24 


0.46 


0.55 


0.51 


0.53 


0.20 


0.20 


0.15 


-0.11 


-0.03 


0.43 


29545 


-1.06 


0.47 


-0.02 


0.44 


0.15 


0.36 


0.29 


0.29 


0.33 


-0.07 


0.00 


0.20 


29598 


-1.06 


0.39 


0.40 


0.48 


— 


0.50 


0.34 


0.26 


0.42 


-0.04 


0.02 


0.28 


29693 


-1.19 


0.26 


0.50 


0.60 


0.64 


0.55 


0.14 


0.23 


0.28 


-0.17 


0.00 


0.36 


29848 


-1.05 


0.54 


0.09 


0.54 


0.64 


0.48 


0.29 


0.32 


0.34 


-0.01 


-0.01 


0.55 


30209 


-0.99 


0.42 


-0.05 


0.46 


— 


0.42 


0.30 


0.31 


0.46 


0.03 


0.05 


0.31 


30345 


-1.06 


0.32 


0.43 


0.44 


— 


0.53 


0.28 


0.28 


0.52 


-0.05 


0.07 


0.36 


30450 


-1.00 


0.40 


0.40 


0.52 


0.51 


0.49 


0.26 


0.29 


0.32 


0.11 


0.05 


0.45 


30452 


-1.00 


0.40 


0.06 


0.50 


0.53 


0.46 


0.27 


0.29 


0.36 


0.02 


0.00 


0.53 


30549 


-1.09 


0.45 


0.13 


0.52 


— 


0.48 


0.34 


0.31 


0.35 


0.00 


0.06 


0.40 


30598 


-1.07 


0.43 


0.05 


0.54 


0.37 


0.54 


0.25 


0.33 


0.41 


0.01 


0.05 


0.52 


30653 


-1.06 


0.47 


0.15 


0.47 


0.46 


0.51 


0.27 


0.21 


0.32 


-0.06 


-0.01 


0.45 


30675 


-1.07 


0.36 


0.41 


0.55 


0.65 


0.50 


0.36 


0.32 


0.42 


-0.02 


0.07 


0.35 


30711 


-1.01 


0.37 


0.32 


0.47 


0.53 


0.45 


0.22 


0.31 


0.41 


-0.07 


0.02 


0.46 


30719 


-1.19 


0.44 


0.42 


0.50 


0.68 


0.46 


0.33 


0.27 


0.48 


-0.09 


0.02 


0.45 


30751 


-1.09 


0.40 


0.32 


0.58 


0.54 


0.48 


0.28 


0.38 


0.33 


0.10 


0.03 


0.42 


30924 


-1.09 


0.48 


0.20 


0.47 


0.48 


0.50 


0.28 


0.30 


0.41 


-0.05 


0.02 


0.40 


30933 


-1.07 


0.44 


0.44 


0.44 


— 


0.58 


0.30 


0.28 


0.34 


-0.08 


0.02 


0.46 


31015 


-1.07 


0.50 


0.00 


0.51 


0.36 


0.42 


0.33 


0.32 


0.40 


0.06 


-0.01 


0.32 


31306 


-1.11 


— 


0.40 


0.41 


0.53 


0.45 


0.29 


0.33 


0.38 


0.03 


0.04 


0.36 


31376 


-1.00 


0.33 


0.43 


0.46 


0.59 


0.43 


0.29 


0.29 


0.34 


-0.01 


0.03 


0.40 


31532 


-1.03 


0.47 


0.38 


0.44 


0.55 


0.47 


0.26 


0.30 


0.39 


-0.08 


0.00 


0.49 


31665 


-1.07 


0.48 


0.02 


0.49 


0.29 


0.46 


0.31 


0.26 


0.35 


-0.06 


0.01 


0.29 


31803 


-1.13 


0.38 


0.25 


0.44 


0.49 


0.51 


0.32 


0.31 


0.28 


-0.10 


-0.01 


0.29 


31845 


-1.06 


0.27 


0.38 


0.41 


0.67 


0.48 


0.27 


0.27 


0.35 


-0.09 


0.02 


0.47 


32055 


-1.12 


0.31 


0.40 


0.55 


0.59 


0.52 


0.27 


0.39 


0.28 


0.03 


0.00 


0.45 


32121 


-0.93 


— 


0.35 


0.44 


0.47 


0.56 


0.29 


0.32 


0.37 


-0.08 


0.02 


0.38 


32151 


-1.07 


— 


0.39 


0.42 


— 


0.49 


0.29 


0.26 


0.34 


-0.04 


-0.01 


0.33 


32317 


-1.07 


0.37 


0.37 


0.54 


— 


0.43 


0.31 


0.33 


0.34 


0.08 


0.03 


0.36 


32347 


-1.09 


0.40 


0.23 


0.51 


0.57 


0.49 


0.27 


0.29 


0.32 


-0.06 


0.03 


0.43 


32583 


-1.08 


0.46 


0.07 


0.46 


0.50 


0.39 


0.31 


0.32 


0.34 


-0.07 


0.00 


0.29 


32627 


-1.10 


0.50 


-0.03 


0.51 


0.36 


0.34 


0.29 


0.26 


0.38 


-0.01 


-0.00 


0.36 


il /UU 


— l.Uz 


U.4y 


U.Uo 


U.ol 


U.jo 


U.4/ 


V.ZZ 


U. J / 


u.jy 


U.Uo 


f\ f\C\ 

U.UU 


U.JZ 


32724 


-1.03 


0.51 


0.11 


0.48 


0.53 


0.39 


0.30 


0.34 


0.39 


-0.04 


0.00 


0.49 


32782 


-1.05 


0.29 


0.46 


0.45 


0.58 


0.41 


0.34 


0.26 


0.27 


-0.07 


-0.03 


0.48 


32871 


-1.01 


0.40 


0.06 


0.62 


0.81 


0.54 


0.29 


0.31 


0.34 


-0.05 


0.01 


0.65 


32874 


-1.14 


0.41 


0.11 


0.55 


0.52 


0.42 


0.26 


0.27 


0.32 


0.02 


-0.02 


0.39 


32933 


-1.13 


0.48 


0.02 


0.56 




0.54 


0.20 


0.29 


0.22 


-0.14 


0.04 


0.53 


32968 


-1.13 


0.41 


0.24 


0.54 


0.50 


0.49 


0.30 


0.25 


0.42 


-0.03 


0.01 


0.48 
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Table 9. continued. 



ID 


[Fel/H] 


[OI/Fe] 


[Nal/FeW 


[MgVFe] NLTE 


[All/Fe] 


[Sil/Fe] 


[Cal/Fe] 


[Til/Fe] 


[TiII/Fe] 


[CrI/Fe] 


[Nil/Fe] 


[BaII/Fe] 


32988 


-1.09 


0.51 


0.02 


0.44 


_ 


0.50 


0.29 


0.29 


0.33 


-0.03 


-0.01 


_ 


33069 


-0.92 


0.39 


0.23 


0.37 


0.59 


0.43 


0.20 


0.34 


0.42 


-0.06 


0.03 


0.42 


33195 


-1.03 


0.48 


0.09 


0.60 


0.66 


0.56 


0.27 


0.32 


0.34 


-0.02 


0.04 


_ 


33414 


-1.05 


0.52 


0.21 


0.49 


0.30 


0.48 


0.32 


0.25 


0.26 


-0.05 


0.00 


0.36 


33617 


-1.09 


0.29 


0.37 


0.54 


0.58 


0.51 


0.32 


0.26 


0.31 


-0.05 


0.02 


0.51 


33629 


-0.98 


_ 


0.31 


0.39 


0.68 


0.55 


0.28 


0.33 


0.37 


-0.02 


0.05 


0.34 


33683 


-1.05 


0.42 


0.00 


0.51 


0.51 


0.42 


0.27 


0.27 


0.39 


0.03 


-0.03 


0.56 


33788 


-1.02 


0.30 


0.37 


0.47 


0.50 


0.43 


0.28 


0.28 


0.33 


-0.02 


0.00 


0.48 


33900 


-1.06 


0.30 


0.47 


0.49 


0.55 


0.46 


0.32 


0.30 


0.39 


-0.05 


0.01 


0.42 


33946 


-1.03 


0.23 


0.32 


0.42 


0.60 


0.45 


0.30 


0.29 


0.33 


-0.03 


0.03 


0.50 


34006 


-1.06 


0.25 


0.44 


0.52 


0.67 


0.48 


0.26 


0.37 


0.27 


-0.08 


0.04 


0.42 


34130 


-1.09 


0.33 


0.43 


0.53 


0.59 


0.51 


0.29 


0.26 


0.34 


-0.04 


-0.01 


0.47 


34167 


-1.10 





0.13 


0.56 


0.48 


0.35 


0.31 


0.25 


0.23 


0.01 


-0.02 


0.04 


34240 


-1.10 


0.47 


0.08 


0.55 


0.55 


0.49 


0.25 


0.35 


0.30 


0.04 


0.03 


0.47 


34502 


-1.08 


0.35 


0.34 


0.47 


0.49 


0.47 


0.28 


0.33 


0.38 


-0.06 


0.02 


0.34 


34579 


-1.08 


0.25 


0.46 


0.52 


0.61 


0.50 


0.28 


0.38 


0.32 


-0.06 


0.02 


0.46 


34726 


-1.01 


0.37 


0.42 


0.49 


0.60 


0.46 


0.28 


0.31 


0.34 


-0.06 


0.03 


0.40 


35022 


-1.08 


0.30 


0.41 


0.52 


0.58 


0.50 


0.33 


0.32 


0.22 


-0.05 


0.02 


0.34 


35061 


-0.99 


0.50 


0.06 


0.50 


0.52 


0.40 


0.24 


0.24 


0.33 


-0.06 


-0.01 


0.39 


35455 


-1.06 


0.45 


0.01 


0.49 


0.49 


0.44 


0.24 


0.23 


0.37 


-0.01 


-0.02 


0.53 


35487 


-1.00 


0.36 


0.30 


0.42 


0.55 


0.44 


0.27 


0.32 


0.43 


-0.05 


-0.01 


0.51 


35508 


-1.05 


0.25 


0.38 


0.41 


0.45 


0.45 


0.27 


0.22 


0.27 


-0.11 


-0.05 


0.46 


35571 


-0.99 


0.45 


-0.02 


0.46 


0.52 


0.45 


0.27 


0.24 


0.34 


-0.09 


-0.04 


0.55 


35627 


-1.08 


_ 


0.04 


0.56 


0.51 


0.38 


0.33 


0.26 


0.38 


-0.09 


0.00 


0.29 


35688 


-1.11 


0.21 


0.40 


0.53 


0.59 


0.47 


0.33 


0.23 


0.34 


-0.04 


0.02 


0.40 


35774 


-1.10 


0.43 


0.24 


0.54 


0.51 


0.53 


0.26 


0.36 


0.31 


-0.03 


0.03 


0.43 


36215 


-1.11 


0.48 


0.18 


0.52 


0.53 


0.53 


0.26 


0.38 


0.34 


-0.04 


0.02 


0.46 


36356 


-1.05 


0.30 


0.26 


0.44 


0.45 


0.50 


0.28 


0.27 


0.43 


-0.02 


0.03 


0.51 


36929 


-1.03 


0.34 


0.45 


0.41 


0.73 


0.50 


0.30 


0.28 


0.38 


-0.04 


0.01 


0.46 


36942 


-0.98 


0.55 


0.04 


0.52 


0.37 


0.45 


0.24 


0.27 


0.37 


-0.16 


0.07 


0.45 


37215 


-1.11 


0.45 


0.25 


0.50 


0.62 


0.48 


0.29 


0.27 


0.47 


-0.07 


-0.02 


0.38 


38075 


-1.07 


_ 


0.42 


0.53 


0.71 


0.54 


0.33 


0.27 


0.31 


-0.05 


0.08 


0.52 


38383 


-1.10 


0.39 


0.11 


0.54 


0.39 


0.44 


0.25 


0.27 


0.26 


-0.06 


-0.01 


0.18 


38399 


-1.08 


0.52 


0.19 


0.55 


0.64 


0.50 


0.29 


0.30 


0.32 


-0.03 


0.04 


0.42 


38896 


-1.02 


0.31 


0.43 


0.48 


0.49 


0.47 


0.28 


0.33 


0.39 


-0.04 


0.01 


0.49 


42490 


-1.07 


0.29 


0.38 


0.48 


0.58 


0.50 


0.24 


0.29 


0.33 


-0.09 


0.03 


0.43 


42620 


-1.09 


0.22 


0.48 


0.49 


0.64 


0.52 


0.29 


0.26 


0.33 


-0.05 


0.03 


0.45 


43370 


-1.05 




0.34 


0.40 


0.51 


0.51 


0.29 


0.28 


0.44 


-0.02 


0.02 


0.31 


44243 


-1.04 




0.41 


0.39 


0.47 


0.53 


0.22 


0.25 


0.49 


-0.04 


0.01 


0.43 


44595 


-1.07 


0.20 


0.44 


0.52 


0.61 


0.46 


0.29 


0.31 


0.38 


-0.00 


0.03 


0.34 


44616 


-1.04 


0.28 


0.44 


0.53 


0.65 


0.46 


0.29 


0.25 


0.33 


0.01 


0.02 


0.37 


45163 


-1.10 




0.46 


0.49 


0.57 


0.54 


0.28 


0.21 


0.35 


-0.03 


-0.00 


0.43 


45895 


-1.04 


0.43 


0.10 


0.57 


0.54 


0.46 


0.26 


0.23 


0.32 


-0.07 


0.00 


0.48 


5359 


-1.03 


0.42 


0.13 


0.55 


0.53 


0.50 


0.31 


0.30 


0.39 


-0.03 


0.03 


0.24 



